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ABSTRACT 

We interpret gamma ray bursts as relativistic, electromagnetic explosions. Specifically, we pro- 
pose that they are created when a rotating, relativistic, stellar-mass progenitor loses much of its 
rotational energy in the form of a Poynting flux during an active period lasting ~ 100 s. Initially, 
a non-spherically symmetric, electromagnetically-dominated bubble expands non-relativistically 
inside the star, most rapidly along the rotational axis of the progenitor. After the bubble breaks 
out from the stellar surface and most of the electron-positron pairs annihilate, the bubble expan- 
sion becomes highly relativistic. After the end of the source activity most of the electromagnetic 
energy is concentrated in a thin shell inside the contact discontinuity between the ejecta and the 
shocked circumstellar material. This electromagnetic shell pushes a relativistic blast wave into the 
circumstellar medium. Current-driven instabilities develop in this shell at a radius ~ 3 x 10 16 cm 
and lead to dissipation of magnetic field and acceleration of pairs which are responsible for the 
7-ray burst. At larger radii, the energy contained in the electromagnetic shell is mostly trans- 
ferred to the preceding blast wave. Particles accelerated at the forward shock may combine with 
electromagnetic field from the electromagnetic shell to produce the afterglow emission. 

In this paper, we concentrate on the dynamics of electromagnetic explosions. We describe 
the principles that control how energy is released by the central compact object and interpret 
the expanding electromagnetic bubble as an electrical circuit. We analyze the electrodynamical 
properties of the bubble and the shell, paying special attention to the energetics and causal 
behavior. We discuss the implication of the model for the afterglow dynamics and briefly discuss 
observational ramifications of this model of 7-ray bursts. 

Subject headings: gamma-rays: burster - magnetic fields 



1. Introduction 

In recent years, a "fireball/internal shock" model of "long" gamma-ray bursts (henceforth GRBs) has 
been developed (e.g. Meszaros 2002; Piran 1999, and references therein). 3 This associates GRBs with 
black hole or neutron star formation during the explosion of rapidly rotating, evolved, massive stars - the 



2 lyutikov@physics.mcgill.ca 
4 rdb3@stanford.edu 

^External shock (e.g. Dermer 2002) and cannonball (Dar &De Rujula 2003) are other proposed models. 



- 2 - 



"collapsar" model. It is proposed that ultra-relativistic jets are formed within the spinning star and that 
these jets are subsequently responsible for the 7-ray emission and the afterglow Woosley et al. (2003). 
This model has been supported by the discovery that GRBs occur preferentially in star-forming regions 
in cosmologically distant galaxies (e.g. Bloom, Kulkarni & Djorgovski 2001), that achromatic breaks (e.g. 
Harrison et al. 2001), indicative of beaming, have been observed in some afterglows and the observation of 
additional luminous components to the late afterglow (e.g. Bloom et al. 2002), which have recently been 
shown to have a supernova spectrum (e.g. Hjorth et al. 2003) in the case of SN 2003dh - GRB 030329. 
The principal phases in this model comprise: 

I Energy Release A source of power associated with a relativistic stellar mass object, now thought to 
be embedded within a star in the case of the long bursts, with luminosity L = 10 50 L S!50 erg s _1 operates 
for a time t s = 100t Si2 s within a region with radius r s = 10 6 r s >6 cm. The energy stored in a combined 
rotational, gravitational and internal form is at least ~ 10 52 L5 i s ,2 erg. The energy release mechanism may 
involve the release of magnetic energy by a torus (e.g. Woosley 1993; Vietri 1998), a nascent magnetar with 
initial angular velocity ~ 10 4 rad s _1 (e.g. Usov 1992; Duncan & Thompson 1992; Thompson 1994; Usov 
1994) or a black hole (e.g. Paczyhski 1986). The magnetic field itself may be produced by strong dynamo 
activity (e.g. Thompson & Murray 2001) or shear (Kluzniak & Ruderman 1998). Interpreting r s as the 
characteristic size of the light cylinder, the associated strength of the magnetic field is B s <~ 10 14 L 1 ^ rj 1 G. 
In an alternative class of models, the energy release involve the formation of a pair plasma by J> 3 MeV 
neutrinos (e.g. Eichlcr et al. 1989). Independent of the source, it is generally supposed that a high effective 
temperature T s = T s , McV - (L / Aureus b) 1,A ~ L^r^ 2 MeV and entropy per baryon (S = 10 6 S 6 k), 
optically thick (tt ~ 10 14 LgQ 4 r~g ) fireball is produced (Cavallo & Rees 1978), varying on a timescale 
ivar ~ r s /c <~ 100r s ,@ (is. Long bursts are argued to have L50 ~ i s ,2 ~ ?"s,6 ~ Sq <~ 1 (e.g. Frail et al. 2001). 

II Flow Formation As the radiation-dominated fireball expands due to "lepto-photonic" pressure, the 
flow is collimated by the surrounding stellar envelope into two anti-parallel jets with opening angle 6 = 
O.10_i. During the subsequent expansion, the energy is converted into ion bulk motion, and becomes 
matter-dominated at a radius, r mat ~ 10 10 cm, where the fluid Lorentz factor saturates with a value Tq = 
10 3 r 0j 3 ~ S(kT s /m p c 2 ) r~j 10 3 ; beyond this radius, most of the energy resides in the kinetic energy of the 
protons. (The stellar photosphere is thought to have a similar radius.) The photons decouple from the 
plasma at a photospheric radius, r p hot which is also in the vicinity of r mat for the envisaged conditions. 

III 7-ray Burst Much of the jet power is dissipated through a series of internal shocks at a radius 
r 7 ~ Tpctvar ~3x 10 12 r§ 3 cm. These shocks are responsible for the re-acceleration of relativistic electrons 
and the production of magnetic field and Dopplcr-shiftcd, 7-ray synchrotron emission, up to <~ GeV energies, 
that is sufficiently well-collimated by the relativistic outflow to escape pair production. This is the Gamma- 
Ray Burst (GRB). The constraint that the highest energy 7-rays be able to escape without producing 
electron-positron pairs, implies that T , 3 ~ 0.3 (e.g. Lithwick & Sari 2001). 

IV Afterglow When r > ct s /9 <~ lO 13 t Si -40Z 1 , the debris takes the from of a shell of cold protons driving 
a blast wave into the surrounding medium with density n = 0.1n_icm~ 3 . An external shock forms and 
when r ~ r rs h ~ 10 16 L^ 2 OZ^^n^i cm, a reverse shock will also form in the exploding debris. The debris 
decelerates after r > rf ree ~6x 10 16 -L5Q 2 tg ^nlj^^ljr^ 1 cm The bipolar blast wave, formed by the shocked 
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circumstellar medium, which now carries most of the energy of the explosion, will further decelerate according 
to T ~ r (r/rf roc ) _3/ ' 2 until it becomes non-relativistic at a radius r nr ~ T^ 3 rf rcc ~ 10 18 cm. During this 
phase, electrons are accelerated and magnetic field is amplified at the outer shock, leading to the formation 
of the afterglow. The non-relativistic blast wave gradually becomes more spherical and evolves to resemble 
a normal, supernova remnant. 

2. Some Problems with the Fireball Model 

The basic fireball model, which we have just sketched, along with its many variations, raises several, 
important questions. Included among these, in temporal order of the flow evolution, are: 

1. How is the entropy of the fireball created? In most models, the release of energy is mediated by 
a strong electromagnetic field which is invoked to create turbulence in an accretion disk (e.g. MacFadyen & 
Woosley 1999), extract the rotational energy of the central black hole (e.g. Kim et al. 2002) and collimate 
and confine the jets (e.g. MacFadyen & Woosley 1999). As the energy release phase lasts for ~ ct s /r s <~ 10 6 
source dynamical times 4 , the magnetic flux must presumably be tied to or trapped by a large, conducting 
mass and the transients should die away so that a quasi-steady, electromagnetic energy flow will be produced, 
similar to what happens with pulsars. The problem is that, if the burst is powered electromagnetically, how 
is the large entropy of a fireball created? As we discuss further below, there is no natural way to accomplish 
this in the vicinity of the source, although there have been suggestions invoking magnetic reconnection in 
an outflowing wind. (Of course, this is not a concern for those models (e.g. Salmonson et al. 2001) where 
the energy release is mediated by neutrinos.) 

2. How are the hypersonic jets made? The fireball model requires that two jets are formed with Mach 
number M <~ 2 1 / 2 To ~ 400 and a ratio of bulk kinetic energy to internal energy ;> M 2 ~ 10 5 . Numerical 
simulations (e.g. Aloy et al. 2002; MacFadyen et al. 2003) and experience with wind tunnels strongly 
suggest that instability and entrainmcnt prevent this from happening inside the star. It is more reasonable 
to suppose that the outflow emerges from the stellar surface (radius i?* = 10 10 i? >ti io cm) with a modest 
Lorentz factor r„, collimated within a cone with opening angle 8 ~ T^ 1 . It will then accelerate linearly 
r cx r due to radiative pressure until either the momentum flux of the radiation field falls below that of the 
ions, or optical depth to Thomson scattering falls below unity. For ion jet the raping radius is (Eq. 22) 
r(r = 1) = 3 x 10 1 2cmL 5 or^jAfi_ 2 - Thus, if an ion jet starts at i?* ~ 10 10 cm with ~ a few, it barely 
has enough optical depth to accelerate to the required T <~ 300. (Note that in this case most acceleration 
happens beyond the photosphere (for a ion jet) at r p h ~ 10 11 cm, see Eq. (15).) 

3. How can the outflow develop large, parallel, proper velocity gradients and avoid producing 
converging streams? 0(1) gradients in the jet proper velocity are invoked in the fireball model in order 
to produce internal shocks, and supply the free energy for particle acceleration. The model is essentially one 
dimensional. It is usually supposed that this variation has its origin in the source which implies that the 
GRB be produced within a radius <; r 2 ci m i n where i m i n is the minimum variation timescale. However, it is 
also envisaged that the jet be collimated and develop ~ (J~0) 2 causally disconnected streams. It is difficult to 
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understand how this collimation can be achieved without producing angular deflections <; T _1 which would 
lead to pair formation by the escaping, high energy 7-rays. (This problem is analogous to the one addressed 
in contemporary cosmology by the theory of inflation with the important difference that, here, it must be 
solved in a continuous flow with a bounding surface.) 

4. What determines the baryon loading of the flow? Independent of the problem addressed in (2), 
in the fireball model, the baryon fraction must be fine-tuned to allow baryons to assume most of the energy 
of the outflow and to attain the large outflow Lorentz factors that are necessary; too small a fraction and 
the energy will escape before the baryon acceleration is complete, too large a fraction and the asymptotic 
Lorentz factor will be too low to allow the highest energy 7-rays to escape (e.g. Meszaros 2002). 

5. Where are the thermal precursors? For ion jets, the escape of thermal radiation at the photosphere 
should produce a thermal precursor with luminosity similar to the main burst (Lyutikov & Usov 2000; 
Meszaros & Rees 2000). These are rarely seen at the 1% level (Daigne & Mochkovitch 2002; Frontera et al. 
2001; Ghirlanda et al. 2003). 

6. How are particles accelerated at relativistic shock fronts? Diffusive shock acceleration, which 
appears to operate efficiently at non-relativistic shock fronts, fails at relativistic shocks, both in the 7-ray 
emitting region and at the external shock front, because only a minority of the back-scattered particles 
can catch up with the advancing shock front. There are promising, kinematic proposals [e.g. Achterberg 
et al. 2001) for addressing relativistic shock acceleration. However, they pre-suppose the existence of a 
subshock in the background thermal plasma and it is not clear how this can be maintained. (Actually, 
thermal and nonthermal particles are not really distinguished in emission models as it is generally assumed 
that a single, truncated, power-law distribution function is transmitted (e.g. Blandford & McKee 1977).) 
More fundamentally, it is by no means certain that relativistic shock discontinuities form at all. It may 
happen that the sharing of momentum between a relativistic outflow and the circumstellar medium happens 
gradually rather than abruptly (e.g. Usov 1994). 

7. How is the magnetic field amplified? In order to produce a high radiative efficiency and fit the 
afterglow light curves, it is necessary for the post-shock magnetic field strength be amplified by a large factor 
over the value it would have due to simple compression. It has been proposed that this amplification is due 
to the Weibel instability (e.g. Medvedev & Loeb 1999). Recent simulations have convincingly shown that a 
long coherence range (much larger than the ion skin depth) of the magnetic field field fluctuation is indeed 
reached Nishikawa et al. (2003); Frederikscn et al. (2003). However, the average values of the magnetic 
energy density, <~ 10~ 3 of the total energy density, is often too low to account for observed synchrotron 
emission. 

8. How is a large degree of 7-ray polarization created? A very high linear polarization (nominally 
80 percent) has been reported in RHESSI observations of GRB021206 (Coburn & Boggs 2003). The obser- 
vation, if typical, is inconsistent with the internal shock model (Lyutikov 2003b) In order to reproduce high 
polarization the internal shock model should make a number of high unlikely assumptions, some of which 
contradict the very fundamentals of the model. There are four assumptions that are made, (i) the field is 
confined to two dimensional plane, presumably the plane of the shock. Magnetic field amplification due to 
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Weibel instability at the shock indeed produces two dimensional fields Medvedev & Loeb (1999); Nishikawa 
et al. (2003); Frederikscn et al. (2003), but the typical size of resulting magnetic structures with linearly 
directed currents is still microscopic, probably tens or hundreds of ion skin depths (which is of the order of 
meters when the fireball is ~ 10 12 cm in size). On larger scales, magnetic field is likely to be three dimen- 
sionally random. In addition, the postshock material must be turbulent: in the fireball model turbulence is 
needed in order to accelerate particle. In order to account for large energy fraction in accelerated electrons 
the turbulent motions should have energy density comparable to the total energy in the shock and thus 
much larger than the energy density in the magnetic field, typically eb < 10~ 3 . This turbulence will easily 
destroy any finely-tuned current structures, (ii) the plane of the turbulent magnetic field is viewed edge-on 
in the rest frame (this requires viewing angle ~ 1/r in the observer frame); (iii) the emitting surface should 
be quasi-planar; this requires that the angular size of the emitting region be A6 < 1/r. (iv) all emitting 
shells must have the same Lorentz factor to be seen edge on (the burst GRB030329 was multi-peaked). As- 
sumptions (ii), (iii) and (iv) are at variance with the fundamental assumption of the internal shock/fireball 
model that every peak is interpreted as being due to collisions of shells with a range of Lorentz factors. 

9. What determines the jet opening angle and its structure? A number of phenomenological jet 
structures have been proposed (e.g. structured, constant or patchy jet). The fireball model neither gives a 
prediction or expresses a preference for the jet structure. 

10. What is the relation between GRBs and X-ray flashes (XRFs)? In the fireball model there 
is no clear relation between GRBs and XRFs, which can be either dirty fireballs, less energetic fireballs, or 
explosions seen " from the side" . 

11. Where are "orphan" afterglows? If the GRB and afterglow emission is associated with jets, as 
described above, then in the most simplistic interpretation, there will be ~ 6~ 2 <~ 1006C 2 "orphan" afterglows 
expected per GRB. Although the current observational constraints are surprisingly poor and the expected 
number is quite model-dependent, it is surprising that no convincing examples have been found so far (e.g. 
Levinson et al. 2002). 

3. Electromagnetic model 
3.1. Overview 

In an attempt to retain the merits of the standard model while addressing these questions, we present an 
alternative, electromagnetic interpretation of GRBs that builds upon earlier models of electromagnetic and 
magnetohydrodynamic explosions (e.g. Blandford & Rees 1972; Benford 1978; Usov 1992, 1994; Meszaros 
& Rees 1992; Ferrari 1998; Kluzniak & Ruderman 1998; Vietri 1998; Spruit 1999; Wheeler et al. 2000; 
Lyutikov & Blackman 2000; Vlahakis & Konigl 2001; Blandford 2002; Lyutikov & Blandford 2002). In our 
specific version of the electromagnetic model, we give a quite different interpretation of the same four phases 
of a GRB introduced in Section 1 (see also Fig. 1). At this point we do not have detail answers to all the 
posed questions; in some cases we offer only a plausible explanation. Also, we specialize to the collapsar 
model although the principles that we describe are easily adapted to other source models that may still be 
needed for the majority of GRBs. 
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I Source Formation (Energy Release) The GRB "prime mover" is a rapidly spinning black hole 
orbited by a massive disk that has just been formed inside an imploding star, or, alternatively, a "millisecond 
magnetar". For qualitative estimates we may use "millisecond magnetar" model by Usov (1992) (see also 
Blandford & Rees 1972), though the numbers will be similar for any relativistic stellar mass source rotating 
with near-critical spin frequency <~ 'O^q^J-i kHz and magnetic field of B s <~ lO 14 ^^ 2 G. The total rotational 
energy, E ~ IVL 2 /2 <~ Lt s ~ 10 52 erg is available to power GRB bursts and the magnetic field is strong enough 
for this energy to be released electromagnetically in a time t s ~ 10 — 100 s. 

We suppose that the outflow primarily takes the form of a large scale Poynting flux and that the 
dissipation rate remains low enough that the power continues to be dominated by the electromagnetic 
component rather than the heat of a fireball well out into the emission region, although there is almost 
certainly an initial phase in which the electromagnetic field is accompanied by a dense pair plasma. 

A rapidly spinning magnetar with a complicated field structure will form a relativistic outflow. The 
behavior of such sources remains an unsolved problem, even in the simpler case of pulsar winds. In this 
paper we adopt a simplifying hypothesis, that the field lines quickly re- arrange to become predominantly 
axisymmctric. Thus we hypothesize that the axisymmetric or "DC" component of the electromagnetic field 
dominates the wave or "AC" component which is either dissipated as heat or diminished through non- 
dissipative rearrangement. In this case the electromagnetic source acts primarily as a unipolar inductor and 
drives a large quadrupolar current flow, rather like what happens in the Goldreich & Julian (1969) model of 
an axisymmetric pulsar. 

II Bubble Inflation (Flow Formation) Initially, the source will inflate an electromagnetic bubble inside 
the star. This magnetized cavity is separated from the outside material by the (tangential) contact disconti- 
nuity (CD) containing a surface Chapman-Ferraro current. This current terminates the magnetic field and 
completes the circuit that is driven by the source. On a microphysical level the current is created by the 
particle of the surrounding medium completing half a turn in the magnetic field of the bubble, so that the 
thickness of the current-currying layer is of the order of ion gyro-radius. 5 

As we show below, the electromagnetic field can be treated as a fluid and behaves similarly to a true 
fluid, with the important difference that the rest frame stress tensor is anisotropic. This allows it to self- 
collimate (for reviews of stationary flow see, e.g. Konigl & Pudritz 2000; Sauty et al. 2002; Heyvaerts & 
Norman 2003). 6 The poloidal and toroidal components of magnetic field are comparable in strength at 
the light cylinder, but the toroidal field dominates beyond this. The velocity of expansion of the bubble is 
determined by the pressure balance on the contact discontinuity between magnetic pressure in the bubble 
and the ram pressure of the stellar material 7 . As the magnetic field strength is strongest close to the 
symmetry axis, the bubble will expand fastest along the polar direction. Eventually the bubble will break 
free of its surroundings and forming a "twin exhaust" along which Poynting flux will flow until either the 
central source slows down or the collimating material itself expands which will both occur naturally on the 
timcscalc t s <~ 100 s. 



5 It is expected that the surface current will be unstable (e.g. Smolsky & Usov 1996; Liang et al. 2003), so that in reality 
the motion of particles will be more complicated and the penetration depth will related to the scale of most unstable modes. 

6 Vlahakis & Konigl (2003a,b) also give examples of collimated MHD outflows applicable to GRBs. 

7 More precisely, since the expansion is supersonic, the pressure balance is between magnetic pressure and the pressure of 
the shocked material, which is of the order of the ram pressure at the forward shock. 
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Outside the star, the bubble will expand ultrarelativistically and bi-conically. After it has expanded 
beyond a radius 

r sh ~ ct s ~ 3 x 10 12 t Si2 cm (1) 

the electromagnetic energy will be concentrated within an expanding, electromagnetic shell with thickness 
<~ r s h and with most of the return current completing along its trailing surface. However, the global dynamics 
of this shell and its subsequent expansion are set in place by the electromagnetic conditions at the light 
cylinder and within the collimation region. 

After break out, the interaction of the magnetic shell with the circumstellar medium proceeds in a 
similar way, except now the velocity of expansion is strongly relativistic. The leading surface of the shell is 
separated by a contact discontinuity (which actually becomes a rotational discontinuity if the circumstellar 
medium is magnetized (Lyutikov 2002a)). Outside the CD an ultra- relativistic shock front may form and 
propagate into the surrounding circumstellar medium. The expansion will still be non-spherical. As long as 
the outflow is ultra-relativistic, the motion is virtually ballistic and determined by the balance between the 
magnetic stress at the CD and the ram pressure of the circumstellar medium. 

The angular distribution of magnetic field (and of the Lorentz factor of the expansion) depends on 
the dynamics of the bubble at the non-relativistic stage and the distribution of the source luminosity. The 
simplest case, which we shall analyze in some detail and which captures the essential features of the outflow, is 
that the outgoing current is confined to the poles and the equatorial plane and closes along the surface of the 
bubble. This produces a toroidal magnetic field that varies inversely with cylindrical radius. Accompanying 
this magnetic field will be a poloidal electrical field so that there will be a near radial Poynting flux, that is 
carrying energy away from the source at almost the speed of light. In addition to the outgoing flux, there is 
a much weaker reflected flux that propagates backward into the flow the information about the circumstellar 
medium. The distribution of reflected current is determined by the outgoing current and the boundary 
conditions. 

III Shell Expansion (7-ray Burst) By the time the shell radius expands to 

r G RB ~ ct s T 2 ~ (Lt 2 J(nm p c 2 ) 1/4 ~3x lO 16 !^ 4 ^ V^cm, (2) 

most of the electromagnetic Poynting flux from the source will have caught up with the CD and been reflected 
by it, transferring its momentum to the blast wave. Simultaneously a strong region of magnetic shear is 
likely to develop at the outer part of the CD (Lyutikov 2002a). Both of these effects are likely to lead to the 
rapid development of current instabilities in the shell that will ultimately result in the acceleration of pairs 
and the emission of Doppler-boosted synchrotron emission in the 7-ray band. Although, we defer discussion 
of the microphysics of particle acceleration to Paper II, we note here that expected radius of GRB emission 
is typically some three orders of magnitude larger than in the fireball model. 

IV Blast Wave Propagation (Afterglow) For r >> tgrb, most of the energy of the explosion will 
reside in the blast wave which will eventually settle down to follow a self-similar expansion. (The structure 
of the energetically sub-dominant electromagnetic shell will also become self-similar.) This is the afterglow 
phase when synchrotron and inverse Compton radiation is emitted throughout the electromagnetic spectrum. 
The initially aspheric expansion will give the appearance of a jet with the "achromatic break" occurring when 
the Lorentz factor becomes comparable with the reciprocal of the observer's inclination angle with respect 
to the symmetry axis. When r > tnr ~ (Lt/pc 2 ) 1 / 3 ~ 2 x lO 18 !/^ 3 ^ 3 ?" 1-1 ^ 3 cm, the blast wave become 
non-relativistic and will become more spherically symmetric, while evolving towards a Sedov solution. 
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3.2. Addressing the Problems of Fireball Models 

Before discussing the dynamical aspects of our model in more detail, we return to the problems that we 
identified with the fireball model and outline how they are addressed in the electromagnetic model. 

1. How is the entropy of the fireball created? Under the electromagnetic model, entropy production 
is deferred until late in the evolution of the explosion where it occurs naturally as a consequence of the 
development of various instabilities. This also addresses the "compactness problem" . 

2. How are the hypersonic jets made? As we discuss further below the effective sound speed is that of 
light and so the jets speeds may be formally subsonic. Jet are collimated naturally through magnetic hoop 
stress. Of course inertial and pressure confinement by a surrounding stellar envelope can also be important, 
thought this is not necessary. 

3. How can the outflow develop large, parallel, proper velocity gradients and avoid producing 
converging streams? No strong constraint need be satisfied because the GRB emission arises at a much 
greater radius than in the fireball model. Furthermore, the emitting region is more strongly coupled causally. 

4. What determines the baryon loading of the flow? As the momentum is carried primarily by 
electromagnetic field, the baryon loading can be negligible. However, it clearly cannot be too large. This 
imposes constraints on the amount of initial loading and entrainment within a stellar envelope. In analogy 
with the Sun one may expect that there are two "phases" within a source: an internal matter-dominated 
one in which large currents arc flowing and an external magnetically-dominated (see Fig. 3). If a flow is 
launched from the magnetically-dominated phase the matter loading may be expected to be small (analogous 
to pulsar wind). 

5. Where are the thermal precursors? The intensity of the thermal precursor is set by the degree of 
lepto-photonic loading which can be arbitrarily small (Lyutikov & Usov 2000; Daigne & Mochkovitch 2002). 
A small precursor seen by Frontera et al. (2001) indeed had very small luminosity. 

6. How are particles accelerated at relativistic shock fronts? As we discuss further in Paper II, 
the particle acceleration for the GRB does not take place at a relativistic shock front but is instead due 
to magnetic field dissipation in the emission region. Electromagnetic energy is "high quality": it can be 
effectively converted into high frequency electromagnetic radiation. For example, in case of Solar flares, the 
primary energy output is non-thermal electrons Benz et al. (2003). The particle acceleration that leads to 
afterglows may be shock-related but could also be due to relativistic MHD modes. 

7. How is the magnetic field amplified? The electromagnetic field is already present and provides 
the dominant energy density during GRB emission. In addition, we suppose that, during the afterglow, the 
magnetic field is supplied by the magnetic shell and may be incorporated into the shocked circumstellar 
gas through interchange instabilities (e.g. impulsive Kruskal-Schwarzschild instability appendix D) or due 
to resistive instabilities of the contact surface. 
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8. How is a large degree of 7-ray polarization created? A strong argument in favor of electro- 
magnetic models comes from the recent report of large polarization in RHESSI observations of the prompt 
7-ray emission from one GRB (Coburn & Boggs 2003). If high polarization is substantiated and found to be 
generic, it would imply that the magnetic field coherence scale is larger than the size of the visible emitting 
region, <~ r/T. Such fields cannot be generated in a causally-disconnected, hydrodynamically-dominated 
outflow. Thus, the large scale magnetic field should be present in the outflow from the beginning and is 
likely to be the driving mechanism of the explosion (Lyutikov et al. 2003). 

9. What determines the jet opening angle and its structure? GRB outflows have large opening 
angles, but do not have a jet in a proper sense. Outflows are non-isotropic so an achromatic break is inferred 
when the viewing angle is ~ 1/r. The jet internal structure corresponds to a " structured jet" with Lg ~ 6~ 2 . 

10. What is the relation between GRBs and X-ray flashes (XRFs)? GRBs are seen from essentially 
all directions in the electromagnetic model. XRFs are GRBs seen " from the side" . The typical total energy 
(inferred from observations of early afterglows) should be similar to GRBs (within an order of magnitude). 
In the only XRF with a redshift this is indeed the case (Soderberg et al. 2003). In a flux- or fluence- limited 
survey, the bursts viewed at large observer inclination should be systematically closer. 

11. Where are the "orphan" afterglows? Since in the electromagnetic model GRB outflows have 
large opening angles the incidence of orphan afterglows should be much less than in the fireball model. 

There are other appealing features of the electromagnetic model, (i) The sources that are invoked are 
very similar to known sources. Pulsar wind nebulae, magnctars and extragalactic jets (as well as, perhaps, 
Galactic jets) are explained as low entropy outflows associated with spinning, magnetized, neutron stars, 
black holes and relativistic disks. What is novel about the GRB is the combination of high field and spin in a 
stellar object, (ii) Variability may be due to the statistical properties of dissipation and not the central source 
activity. Magnetic fields arc non-linear dissipative dynamical system which often show bursty behavior with 
power law PDS. (iii) "Standard candle" - the narrow distribution of GRB energies (inferred from prompt 
emission, Frail et al. (2001), from afterglows, Panaitescu & Kumar (2002), and from K a lines, Lazzati 
et al. (2002)) may be related to the total rotational energy of a critically rotating relativistic object - a 
one parameter (mass) family, (iv) Correlations of GRB properties: hard-to-soft temporal evolution of GRB 
spectra and -Epeak oc \fL correlation naturally occur in the model (see section 11). 



The electromagnetic model exhibits some simple generalities that derive from treating it as a circuit. 
Suppose that the source is threaded by a magnetic flux, $ ~ £? s r 2 <~ 10 26 G cm 2 , and has a typical angular 
velocity n ~ 10 4 rad/s" 1 . The source will generate an EMF, V, and an associated current, / 



3.3. GRBs as electromagnetic circuits 





There will be an associated current 
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where B s = 10 14 £> S 44 G is the source magnetic flux density and ^i oa d ~ 100 is the total impedance of the 
source and the emission region, which is of order the impedance of free space under general electromagnetic 
and relativistic conditions. The source region can be thought of as a generator capable of sustaining an EMF 
£. (Under most conditions the maximum energy to which a particle of charge Z can be accelerated will be 
limited by ~ Ze£ ~ 3 x 10 22 eV, way above the highest energy of the observed cosmic rays of 3 x 10 20 eV.) 
This implies that the power dissipated in the load is L ~ £ 2 /Zioad- The load consists of external medium, 
against which PdV work is done by the expanding shell, and radiation (some of the energy of the shell is 
radiated away as a prompt emission). An equivalent and useful way to think about this is to say that there 
is a strong, quadrupolar current distribution outward along the axes and inward along the equator (or vice 
versa). Our proposal differs from the conventional interpretation principally through the assumption that 
the current flows all the way out to the expanding blast wave, rather than completes close to the source (c.f. 
Fig. 2). 

We now consider, in more detail, the dynamics of the expansion. We divide the explosion into four 
phases, source formation (t << t s ), bubble inflation (r s /c << t <, t s ), shell expansion (t s <,t <J tgrb/c and 
blast wave propagation (rcRB/c ^5 t). Each of these phases involves distinct, dynamical behavior requiring 
different approximations to describe. 

4. Source Formation 

4.1. Nature of the Compact Object 

As explained above, there is now good circumstantial evidence linking at least some GRBs with simul- 
taneous supernova explosions - the collapsar model. In its original and most common form, (e.g. Woosley 
1993), it is supposed that the core collapse of a massive star leads to the formation of a massive black hole 
orbited by a dense, thick accretion disk 8 It is also assumed that magnetic flux is generated locally by dynamo 
action (e.g. Thompson & Murray 2001). The actual power may derive from the spinning spacetime of the 
black hole or at the expense of the binding energy of the orbiting gas. For a rapidly spinning hole of mass 
<~ 10 M Q , a field of strength B s <~ 10 14 Ll' o 2 G suffices in either case. The simplest descriptions of these pro- 
cesses comprise solutions of Maxwell's equations in a curved spacetime under the force-free approximation. 
They describe, at least conceptually, a stationary, axisymmetric flow of electromagnetic energy and angular 
momentum away from the surface of the hole and the disk. There is an associated current distribution that 
is quadrupolar so that the sign of the radial component of the current changes with latitude. 

Real source is unlikely to be either stationary or axisymmetric. The field configuration may well be 
unstable and irregularities in the disk will break the symmetry. Our fundamental assumption, that underlies 
most of what follows, is that these instabilities do not develop to large nonlinear amplitude and completely 
disrupt the outflow. In other words, the large scale field is simple and approximately axisymmetric beyond 
the light cylinder. This large scale regularity can come about either because the smaller scale magnetic 
structure is erased through dissipation or, in the case of a central neutron star, through a topology-preserving 
re-arrangement of the magnetic flux. Our "DC" model is therefore rather different from other "AC" proposals 
that GRBs be powered by a rapidly varying Poynting flux (e.g. Lyutikov & Blackman 2000; Spruit et al. 
2000; Sikora et al. 2003). 



The observed coincidence between the burst and supernova is contrary to the expectation of the alternative supranova 
model (Vietri 1998). 
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In the magnetar model, it is proposed that the neutron star is born with a rotation frequency ~ 
3LgQ 2 ^2 2 kHz and a dipole field of strength B s <~ 10 14 L 50 1//2 t~ 2 1 G so that it can produce the inferred, 
electromagnetic power and energy. If the dipole is inclined with respect to the rotation axis, we expect that 
some of the open magnetic flux from the northern magnetic pole finds its way into the southern hemisphere 
beyond the light cylinder and vice versa. Adopting our conjecture, provided that the dipole is not too 
inclined, the asymptotic electromagnetic configuration is roughly axisymmetric with field lines ending up in 
the hemisphere from which they started. In other words, the corrugations in the current sheet that separates 
the north seeking and south-seeking field lines smooth out, with little dissipation, close to the light cylinder. 
9 Much of what follows is predicated on this hypothesis, that a large scale quadrupolar current flow is 
established and provides a good description of the subsequent evolution of the bubble and the shell {e.g. 
Blandford 2002). 

In the immediate vicinity of the source the plasma is separated into two phases: matter-dominated and 
magnetically dominated (c.f., Solar photosphere, pulsar magnetospheres). Superstrong magnetic fields are 
generated in the dense medium (a disk or a differentially rotating neutron star-like object). Buoyant magnetic 
field lines emerge into the tenuous magnetosphere while remaining anchored in the matter-dominated phase. 
Strongly relativistic outflow is generated in the magnetically dominated phase. 

Although the magnetosphere is comparatively small, the electrodynamical conditions at the light cylin- 
der constitute a boundary condition for the eventual, relativistic outflow, much like what happens with 
vacuum electromagnetic radiation. We can think of these conditions either as establishing the current dis- 
tribution or, equivalcntly, as defining the subsequent evolution of the electromagnetic field. As the source 
will, typically, remain active for <~ million dynamical times, we suppose that it will be able to settle down 
quickly to a quasi-steady state evolving slowly as the hole or neutron star slows down on a time scale of <~ 
100 s. 



4.2. Dissipation at the Source, r <~ r s 

In any scenario some fraction of the central source luminosity is likely to be dissipated close to the 
source. In other words, there is a source impedance Z s . In the case of electromagnetic extraction of energy 
from a spinning black hole, this dissipation occurs beyond the event horizon. For a magnetar, the neutron 
star impedance is negligible and Z s is dominated by what happens in the magnetosphere. In the fireball 
model, it is implicitly assumed that all of the energy released is quickly converted into heat, forming a high 
entropy per baryon, thermal plasma. In other words, the load impedance Z^ is located close to the source. 
By contrast, in the electromagnetic model this does not happen and the energy flows way from the light 
cylinder mainly in the form of an electromagnetic Poynting flux and the load impedance is located in the 
emission region. We argue that this is likely to be the case because, somewhat paradoxically, it becomes 
harder to convert electromagnetic energy directly to pair plasma the stronger the magnetic field becomes. 
The reason is that the plasma surrounding the source is highly conductive. In such plasma there is usually 
plenty of charge available to screen the component of the electric field along the magnetic fields, so that the 
first electromagnetic invariant is close to zero: E • B ~ (the perpendicular component of the electric field 
just defines a plasma velocity as long as E < B, see below). 



9 There is a good precedent for this behavior in Ulysses observations of the quiet solar wind (McComas et al. 2000) which 
reveal that, despite the complexity of the measured surface magnetic field, the field in the solar wind quickly rearranges to form 
a good approximation to a Parker (1960) spiral. 
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In strongly relativistic plasmas, a possible source of (inertial) resistivity is related to the break down of 
the E • B <~ approximation in cases when the real charge density falls below some critical value. In the 
case of rotating magnetic field, this minimum charge density needed to short out the component of electric 
field along the magnetic field is known as the "Goldreich- Julian" density 

n G J = -z 10 17 cm- 3 5 

Z7rec 

The minimum energy density associated with this amount of plasma is najmc 2 . It is convenient to introduce 
a parameter a as the ratio of the magnetic energy density tig = B 2 /87r to the total plasma energy density 
u p including photons (u p has contributions from rest mass of ions and pairs plus internal energy) 

If the energy density u p is dominated by leptons with density given by Eq. (5), then the ratio (6) is roughly 
the ratio of the light cylinder radius to the electron gyro radius which can be as large as <~ 10 18 . 

A conservative flow of plasma may not be able to satisfy the constraint that the local charge density 
always exceeds uqj. If this happens, a gap will develop where the field-aligned electric fields are nonzero 
E • B =^ 0. However, the maximum potential drop that is available for dissipation will be limited by 
various mechanisms of pair production. Typically, after an electron has passed through a potential difference 
AV ~ 10 9 — 10 12 V it will produce an electron-positron pair either through the emission of curvature photon 
or via inverse Compton scattering. This will be followed by an electromagnetic cascade and the newly born 
pairs will create a charge density that would shut-off the accelerating electric field. The typical potential 
difference required for the creation is orders of magnitude smaller than the total available EMF £ ~ 10 22 V. 
In this case the magnetization parameter may be estimated as 

-^~10 10 (7) 

Another possible way through which an electromagnetically-dominatcd flow can create entropy directly 
and reduce a appreciably is through the development of an electromagnetic turbulent cascade operating down 
to wavelengths, A m ; n , so small that electromagnetic energy can be dissipated directly in particle acceleration. 
This is analogous to the viscous dissipation that terminates a fluid turbulence spectrum. If this really can 
operate then it is hard to see how more than a few percent of the electromagnetic energy will be dissipated 
in this fashion. More specifically, for a differential cascade a typical cascade time is the interaction time 
multiplied by the logarithm of the outer and inner scale (Zakharov et al. 1992). If the typical interaction 
time is a light crossing time, then we can estimate a ~ ln(r s /A m ; n ) ~ 100 (where A m i n of the order of Larmor 
radius), so that the flow will remain electromagnetically dominated. 

Another source of dissipation is magnetic rcconncction. This is surely important if the outflow retains an 
AC component, contrary to our hypothesis, or if current sheets develop as a result of the nonlinear evolution 
of electromagnetic instabilities. Driven relativistic reconnection may proceed at velocities approaching the 
speed of light (Lyutikov & Uzdensky 2003), however the initial development of dissipative current sheets 
occurs on a timescale intermediate between very long resistive and very short dynamic (light travel) timescale 
and may be too slow to dissipate much of the magnetic energy density (Lyutikov 2003a). 

Thus, an abundance of pairs can be created without changing the electromagnetic dominance and the 
more pairs are present, the harder it becomes to create the small "gaps" that would be necessary to replenish 
them. Put another way, the pair density required to supply the electrical current and space charge scales 
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linearly with the field strength, while the electromagnetic energy density scales as its square. The stronger 
the field, the more likely it is to persist into the outflow. It is because GRBs are so powerful that the 
dissipation in the source is probably low. 

5. Bubble Inflation, r < ct s 

In the previous section we have argued that the flow formation occurs on a scale of a light cylinder, 
r ~ r s <~ 10 6 cm and that the dissipative processes are not likely to drain all the potential EMF, so the 
flow is likely to remain magnetically-dominated. The velocity with which the plasma leave the neighborhood 
of the light cylinder is determined by the details of the acceleration and magnetic field structure at the 
source. For the purpose of this paper we leave the question of the detail structure of the central source open, 
assuming that the formation of the flow occurs in a way similar to pulsar magnetosphcres (Michel 1969; 
Goldreich & Julian 1970), with an important difference - in GRBs pressure effects may play an important 
role. Beyond the light cylinder the magnetized flow generated by the central source will expand due to 
magnetic and pressure forces (see also Wheeler et al. 2000; Moiseenko et al. 2003). It will become super- 
fast-magnetosonic and thus causally disconnected from the source (Goldreich & Julian 1970). For subsequent 
evolution the conditions close to source may be considered as boundary conditions at which the rate of energy 
and magnetic flux injection is some given function. 

The asymptotic structure (at r » r s ) of axisymmetric, magnetized outflows is a challenging problem 
that has a considerable literature (e.g. Heyvaerts & Norman 2003, and references therein). Heyvaerts and 
Norman argue that a perfect, non-relativistic MHD flow with five conserved constants of the motion evolves 
either to a state where the current is confined to the axis and a finite number of thin current sheets or 
that the current closes and that the outflow energy flux becomes purely mechanical (for relativistic flows the 
latter happens at distances much larger than astrophysically relevant scales). In the context of a quadrupolar 
current distribution, the first option is equivalent to stating that the current becomes concentrated on the 
axis, surface and equator. They also argue that the second possibility is ultimately favored but that finite 
flows may not, in practice, achieve this state. 10 . 

As the wind expands, it starts to interact with the surrounding medium, so that its properties are 
determined both by initial conditions and interaction with the surrounding medium. The wind is slowed 
clown by interaction with the stellar material, so that initial expansion is non-relativistic, and later, after 
breakout, it becomes strongly relativistic. The physical processes governing these two phases are quite 
distinct. We consider them in turn. 



5.1. Non-relativistic expansion 

Consider a newly-formed, compact object inside a star or other dense gas distribution generating an 
EMF £{t) and driving a quadrupolar current flow I(t) as described above. This will inflate an electromagnetic 
bubble expanding at a rate controlled by the external gas density. Let the bubble radius be R(6, t) where 
9 is the polar angle measured from the symmetry axis defined by the spin of the compact object. For the 
moment, suppose that the bubble expands non-relativistically. 



10 An alternative way to justify this particular electromagnetic configuration is to note that it minimizes the magnetic energy 
associated with a given quantity of magnetic flux. 
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We expect magnetic flux (integrated over the meridional plane) to cross the light cylinder and be 
supplied to the bubble at a rate 4> ~ /j,oIc/2tt. Similarly, Poynting flux will be supplied to the bubble at a 
rate U EM ~£I U . We can also compute the magnetic flux <f> = CI and the energy stored within the bubble 
Uem = CI 2 /2, using the self inductance C. If, as we discuss further below, the magnetic field in the bubble 
is predominantly toroidal between cylindrical radii n7 m i n and tu max = Rsin9, this is given by 

fdzln(^). (8) 
2n J \ tu min / 

We therefore see that, if the bubble expands homologously and sub-relativistically, the rate of supply 
of both flux and energy exceeds the rate at which the flux and energy can be stored by a factor <~ 
[ln(ra7 max /tn m i n )((iz/(it)/c] _1 (c/. Rees & Gunn 1974). Therefore, too much flux and energy is generated by 
the source when the expansion speed is less than ~ 0.1c. 

The fate of this surplus flux and energy depends upon the amount of dissipation within the bubble. If, 
despite strong driving, the resistance in the electrical circuit is sufficiently low (much less than ~ 100f2), then 
the electromagnetic energy would be either reflected back to the source changing its properties (this can and 
does happen with solid conductors, e.g. waveguides, but is unlikely to happen in a plasma environment), or 
there will be "inductance breakdown", so that C will decrease below the estimate (8). This will be achieved 
by destroying axial symmetry of the flow by MHD instabilities and creation of smaller scales current, much 
along the lines of what has been suggested in the Crab Nebula by Begelman (1999). On the other hand, if 
the resistance in the electrical circuit is sufficiently high, magnetic flux and electromagnetic energy will flow 
toward those parts of the current flow where the resistance is located. Magnetic energy then can be dissipated 
at a very high rate (this is similar to the case of driven magnetic reconnection, which in relativistic plasma 
may proceed at the speed approaching a speed of light (Lyutikov & Uzdensky 2003)). The total electrical 
resistance needed can be estimated from the ratio of the potential difference at the light cylinder to the 
current as <~ 60f2 times a logarithmic factor and we will estimate the total resistance as lOOfi. In this paper, 
we assume that the current flows out to the boundary of the bubble and the Poynting flux is transformed 
into heat until the expansion speeds exceeds <~ 0.1c. Thereafter, it can be accommodated non-dissipatively 
by the expanding, electromagnetic bubble. 

However, if the resistance in the outer part of the bubble always exceeds <~ lOOfi, then the current will 
complete closer to the compact object and the outer parts of the bubble will comprise hot radiation/pair- 
dominated plasma. This is the implicit assumption underlying fireball models of GRBs that invoke electro- 
magnetic sources. 

As mentioned above, ideal MHD flows are naturally collimating. At the non-relativistic stage of expan- 
sion there is another collimating effect due to resistive dissipation, which, as we have argued, must be present 
in the flow. Most of the dissipation is likely to occur near the axis where the current density is highest and 
the susceptibility to known instability is the greatest. In this case a lateral flow of energy will set in carrying 
the poloidal field lines with it towards the axis (dissipation of magnetic energy on the axis will result in a 
loss of magnetic pressure, which resists the inflow of plasma towards the axis, and will be communicated to 
the bulk of the flow by a rarefaction wave propagating away from the axis). This, in turn, will leads to the 
pile- up of magnetic field near the axis and to faster radial expansion (the toothpaste tube effect). 

Next we consider the structure of the axial current, whose radius tz7 m ; n we have already introduced. If the 
energy that is dissipated were to be radiated immediately and to escape freely, then the field structure would 



If there is outward and return current at intermediate latitude, then E/em = / d£I, where I is the total enclosed current. 
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evolve to become force- free everywhere. This would imply that zu m i n ~ r s , the radius of the light cylinder. 
However, this is surely not the case. The optical depth of the plasma will be far too large and the radiation 
will be trapped and thermalize. This implies that pressure is likely to become quickly important in the core. 
For r » r s , the net poloidal field is quite small. The simplest structure to consider is that of a Bennett 
pinch. In a Bennett pinch, the internal energy associated with the plasma is 'S^qI 2 /16tt per unit length so 
that the mean value of a is 41n(n7 max /n7 m i n )/3. A substantial plasma energy density is needed to oppose the 
magnetic stress independent of the choice of w m i n . However Bennett pinches are notoriously unstable and so 
this is also unlikely to be a complete description of the current. The instabilities that would develop within 
the core lead to a less ordered and dynamic magnetic field with pair plasma contributing to the overall stress 
tensor to an extent controlled by the balance between its rate of creation and annihilation which we discuss 
below. Furthermore, we suppose that these instabilities lead to dissipation and that the effective impedance 
in the circuit automatically adjusts to match that required to account for the electromagnetic energy and 
flux supplied by the source at the light cylinder and the rate of expansion of the bubble. 

In summary, we adopt a particularly simple distribution of currents generated by the central source: 
along the axis, the surface of the magnetic bubble and closing in the equator. More general current distri- 
bution do not change the picture qualitatively, but will lead to quantitative changes in our conclusions. 



5.2. Dynamics of magnetic bubble inside a star 

Electromagnetic bubbles are crucially different from expanding fireballs in another way. They are 
naturally self-collimating and will create bipolar outflows even if the stellar mass distribution is spherically 
symmetric. 

The dynamics of such non-sphcrically expanding bubble may be described using the method of Kom- 
paneets (1960) which was initially developed for propagation of non-spherical shocks (for astrophysical ap- 
plications see Icke 1988; Bisnovatyi-Kogan & Silich 1995). Consider a small section of non-spherical non- 
relativistically expanding CD with radius R(t, 9). The CD expands under the pressure of magnetic field so 
that the normal magnetic stress at the bubble surface is balanced by the ram pressure of the surrounding 
medium. At the spherical polar angle 9 the CD propagates at an angle 

OlnR 

tana " — w~ () 

to the radius vector. Balancing the pressure inside the bubble B 2 /(8ir) — I 2 /(2irc 2 R 2 ) with the pressure of 
the shocked plasma we find 



dRy i 2 (t) 



dt J e 2irR 2 sin 2 9p(R, ( 



1 



d\nR 
89 



(10) 



where k is a coefficient of the order of unity which relates the pressure at the CD to the pressure at the 
forward shock. 

Equation (10) shows that non-spherical expansion inside the star is due both to the anisotropic driving 
by magnetic fields and collimating effects of the stellar material (the term in parenthesis, which under certain 
conditions tends to amplify non-sphericity) . Note that this use of the Kompaneets approximation assumes 
that the shock and the CD are located close enough so that there is no lateral (in 9 direction) redistribution 
of pressure in the shocked material. This is a good approximation for accelerating shocks and is used here 
only to illustrate qualitative behavior of the solutions. 
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The rate of expansion of the bubble inside the star depends upon the dynamics of the stellar envelope 
and the time evolution of the current I(t). For a given dependence p(R,0) and I(t) Eq. (10) determines 
the velocity of the CD. Generally solutions will be strongly elongated along the axis. A simple analytical 
solution for I, p ~ const is 

^M) = (-4V /4 ^- (ID 

\ 7r pc ) sin a 

(current is related to the luminosity by Eq. (4). Qualitatively, the bubble and the forward shock will cross 
the iron core (r c <~ 2.5 x 10 8 cm) in t <~ r cy /p/B(r c ) ~ .3 sec, short enough to produce an ample supply of 
56 Ni Woosley et al. (2003). If we define M(R) as the stellar mass external to radius R, then we can also 
estimate the breakout time 

Mo) (r^J L » 8 (12) 

The electromagnetic bubble can be confined equatorially by the star for the duration of the burst ibreakoutC 7 '"/^) 
100s and will expand non-relativistically as we have assumed. However the expansion along the axis proceeds 
on a short timescale and breakout should occur early in the burst. Furthermore, at the time of breakout, 
the axial expansion speed of the bubble will become rclativistic. 

Non-relativistic expansion lasts for several seconds along the axis. This time is much shorter than 
the burst duration. Recall that when the velocity of expansion is < 0.1, a considerable fraction of the 
magnetic energy is indeed dissipated. But since the flow quickly becomes relativistic, the relative fraction of 
dissipated energy is small, so that the flow magnetization remains large, a > 100. An important advantage 
of the electromagnetic models is that as long as a > 1, the asymptotic evolution of the flow is mostly 
independent of a. Most of the dissipation described above will result in creation of lepto-photonic plasma, 
which decouples after photosphere, so that the remaining flow remains strongly magnetized. For a > 1 the 
energy associated with the thermalized component is not important for the flow dynamics. 

In summary, our contention is that, initially, when the bubble expands non-relativistically, the dissipa- 
tion is concentrated along the pole and the current returns to the source along the surface of the bubble. 
The polar current pinch is supported against collapse at its core by a combination of plasma pressure and, 
possibly, by dynamical, disorganized magnetic field. After breakout, the expansion becomes relativistic and 
the resistance falls so that the electromagnetic energy that is still being supplied by the source is mostly 
absorbed by the inflating bubble and by doing work against the surroundings. 



5.3. Early optically thick expansion: mini-fireball 

In the previous subsection we have argued that energy release and initial non-relativistic stage of ex- 
pansion are necessarily accompanied by partial dissipation of the magnetic energy, but the flow is likely to 
remain magnetically dominated with it> 1. For any reasonable a < 10 10 the lepto-photonic component will 
be optically thick near the central source and consequently in a thermodynamic quasi-equilibrium. In this 
subsection and in appendix A we discuss the optically thick, quasi-spherical expansion of a relativistically 
hot, magnetically-dominated flow after the flow became weakly relativistic and no further dissipation of 
magnetic energy is happening in the flow, but at the early enough stages so that photons remain trapped 
(see also Vlahakis & Konigl 2003a, b; Fendt & Ouyed 2003, for a more extensive analysis). 

Under the electromagnetic hypothesis, most of the energy released by the source comes out in the form 
of Poynting flux. However, as we argued above, there must be some dissipation that would lead to creation 
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of a lepto-photonic component. In addition, some ions may be present in the flow. The luminosity of the 
source then can be written as 

L = j dnr 2 r 2 (3c(b 2 /2 + w) (13) 

where w is plasma enthalpy, b is a toroidal magnetic field in the plasma rest frame, Q is a solid angle. The 
luminosity (13) includes contributions from the rest mass energy density of ions and pairs, their kinetic 
pressure and the trapped photon gas. At early stages the plasma enthalpy is strongly dominated by lepto- 
photonic plasma with a temperature 

T ( L 

~ \aAnr 2 l3cT 2 (l + a) 

where AO is a typical opening solid angle. (The magnetization parameter in this case is <r = b 2 /2/w). 
For any reasonable values of a close to the source this temperature is high enough, so that pairs are freely 
produced. At breakout, rbreakout ~ 10 10 cm, plasma is moving weakly relativistically, F <~ 1, so that the 
temperature may still be high enough for plasma to be pair-dominated 

Tbreakout ~ 100 keVI^O:^ 17 Vr^L«t,io ( 15 ) 




At this point the flow will accelerate to relativistic velocities. Initially, the expansion is mostly pressure- 
driven, even in the strongly magnetized case. This results in dynamics qualitatively similar to the unmag- 
nctized case. During outflow, the wind plasma accelerates r ~ r while its density, pressure and temperature 
decrease n ~ r~ 3 , p <~ r~ 4 , T ~ r _1 . During the pressure-driven expansion the flow becomes superfast 
magnetosonic T 2 > a, while the magnetization parameter remains approximately constant (appendix A). 
(Note that the magnetization parameter is equal to the ratio of the Poynting to the particle fluxes. Since it 
remains approximately constant for quasi-spherical expansion, there is little transfer of energy between the 
magnetic field and the plasma at this stage for the assumed quasi-spherical geometry of the outflow.) 

When the temperature falls below <~ 10 — 20keV, most of the pairs annihilate. This occurs at r p h ~ 
10 11 cm. This suddenly reduces the optical depth to Thomson scattering below unity. (Under certain 
conditions photons may remain trapped in the flow (Section 5.4 below, also Lyutikov & Usov 2000). In 
this case, thermal driving by photon pressure continues, until the thermal photons escape.) As a result the 
lepto-photonic part of the flow decouples from the magnetic field and a increases by roughly seven orders 
of magnitude to a <~ 10 9 . By this time most of the thermal energy has been spent on accelerating the flow 
to T ~ 10. Beyond the photosphere, thermal photons propagate freely. The thermal radiation from the 
lepto-photonic component has a rest-frame temperature T ~ 10 — 20keV times a boost due to the bulk 
motion. This thermal radiation, which should peak around ~ 100 keV may put constraints on the initial 
a (Lyutikov & Usov 2000; Daigne & Mochkovitch 2002). There are indications that the thermal precursor 
has been observed Frontera el al. (2001), with intensity <~ 1% of the total burst intensity. This can be used 
to estimate the magnetization parameter below photosphere as a <~ 100, so that beyond the photosphere 
a ~ 10 9 . 



5.4. Thomson and pair production depths 

Next, we consider conditions on the optical depth to Thomson scattering and pair production in electro- 
magnetic models beyond the photosphere when the ejecta plasma is cold, collisionless (so that the number 
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of pairs is conserved) and strongly magnetized (a ~ 10 9 ). Under certain conditions, determined below, it 
may remain optically thick to Thomson scattering and to pair production. 

Under the approximation above, the total energy is carried by the Poynting and particle fluxes (the 
latter includes a contribution from the pair rest mass and proton rest mass). If we neglect possible effects 
of pair production due to the dissipation of magnetic field, the magnetization parameter becomes 

° 47rc 2 m±n' ± (l + { ' 

where £ = m i n' i /m±n' ± is the ratio of ion to pair mass fluxes. (In a standard fireballs, a = 0, and pair arc 
dynamically unimportant £ <~ m p /m e ; for magnetized, pair-loaded flows , a » 1, £ < 1. ) The total energy 
flux then can be written as 

L = rWc^±^AQ (17) 
a 

and we find that 

n '± = A^r 2 mc 3 r 2 (l + ( r)(l+0 (18) 
Introducing the compactness parameter 

_ L<j T r< 

c Attr 2 r 2 mc 3 { ' 

where r' ~ r/T is the typical size of emission region in the rest frame, we find that the optical depth to 
Thomson scattering is 

TT = ^irl = l s (20) 

T Afir 2 mc 3 r 2 (l + a)(l+0 (1 + <t)(1+0 1 ' 

In order to escape, we require that photons must have Tr <1. This places a lower limit on the radius: 



/ Lot?' 

T> y A^mc 3 r 2 (l + a)(l+0 ( ' 

which under the assumption r' ~ r/T gives 

r >= (22) 

Afimc 3 r 3 (l + ( r)(l+e) 1 j 

Thus, for large magnetization, a > 1, even for pair dominated flow, 1 the ejecta become optically thin 
to Thomson scattering at 

r > 10 8 (1 + a)g 1 AQzlTl cm (23) 

(For an electron-ion plasma £ = m p /m e , so that for a given a the flow becomes optically thin at radii which 
are three orders of magnitude smaller than for a pair plasma.) 

Next, we estimate the radius at which the optical depth to pair production falls below unity. Assume 
that most photons emerge with energy E\, r ~ 200 keV. The optical depth to pair production then becomes 



= n>(E br )^r> = ^§^^ Tf ^ (24 , 



where e 7 is the fraction of the luminosity emitted in 7-rays, (7 7 _ 7 = (11/180)<7t is the cross-section for pair 
production by particles with power law spectrum (Svensson 1997), and primes denote quantities measured 
in the outflowing frame. 
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Using the definition of the compactness parameter (19), we can write the optical depth to pair production 
(24) as 

180 " 7 c E br 

This can be called photon compactness. Since mc 2 / E\, T ~ 2, 



r 7 = T^ e 7 r ^— (25) 



t 7 - 0.1H C (26) 

The condition r 7 _ 7 < 1 gives 

r ->slr^= 3 * 10,6i » r = 2 ™ <27) 

The ratio of t t to r 7 may then be used to find the ratio of photons to electrons 



'7 

n: 



^ ~ 0.1IV(1 + » 1 (28) 
n± 

Thus, there are many more photons than electrons and the flows will first become optically thin to Thomson 
scattering, but will remain optically thick to pair production up to much larger distances. For strong 
magnetization, a ^> 1 the flow will become Thomson thin right after breakout, while remaining optically 
thick to pair production up to r ~ 10 16 cm (Eq. (24)). Note that in the electromagnetic model pair 
production may play a more important role than in the hydrodynamic models since it can regulate the 
component of electric fields along magnetic field and control acceleration of pairs. 



5.5. Magnetic acceleration and collimation in the relativistic regime 

The behavior of magnetized winds depends both on the conditions at the source (e.g. the lateral distri- 
bution of the energy and magnetic fluxes) and on the subsequent dynamics. Since for t < t s the expansion is 
quasi-stationary, we can assume that, at this stage, the progenitor produces a steady relativistic wind. There 
is an extensive literature on acceleration and magnetic collimation of outflows from pulsars (Michel 1971; 
Bcnford 1984; Sulkanen & Lovelace 1990; Begelman & Li 1992; Fendt et al. 1995; Bogovalov & Tsinganos 
1999; Lyubarsky 2002) and from disks (Blandford & Payne 1982; Lovelace et al. 1987; Contopoulos & 
Lovelace 1994; Camenzind 1995; Fendt & Ouyed 2003). (Unfortunately, even in the cleanest case of pulsar 
winds, the dynamics of relativistic MHD outflows remains an unsolved problem). 

Magnetic outflows can be naturally self-collimating through the action of magnetic hoop stresses, oc 
B^/^tt. Generally, this magnetic hoop stress is counterbalanced by the gas pressure and magnetic field 
gradients, so that the flow evolution (collimation or de-collimation) depends on this balance (e.g. Heyvaerts 
& Norman 1989; Chiueh et al. 1991; Vlahakis & Konigl 2003a; Fendt & Ouyed 2003). Collimation also 
affects the magnetization of the flow. For radially expanding flow, a remains constant since both the magnetic 
energy density and the plasma rest mass energy density scale as cx r~ 2 , so there is no transfer of energy 
between them (in this case magnetic hoop stresses are exactly compensated by the gradient of magnetic 
energy density). Collimation (or decollimation) leads to transfer of energy from magnetic field to plasma (or 
vice versa). 

An important property of ultra-relativistic outflows is that their motion is virtually ballistic, so that 
any collimation should be achieved close to the source where the flow is only mildly relativistic. This is a 
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wcll-known problem for stationary pulsar outflows and in AGNs (e.g. Begelman & Li 1992; Bogovalov & 
Tsinganos 1999). The reason for the weak collimation is that far from the source the gas pressure and the 
poloidal magnetic fields are unimportant and the magnetic stress oc is almost exactly balanced by the 
oppositely directed electric stress oc (3 2 B 2 . Thus, the resulting stress is quite small oc B 2 /T 2 (e.g. Bogovalov 
2001). Increasing the magnetic field strength does not help either, since in the strongly magnetized outflow 
the rest frame energy density is oc B 2 , the resulting lateral acceleration is independent of the magnetic field 
strength. A simple way to see this is to consider the MHD force balance equation in the 6 direction (see Eq. 
(C3)) 

d t l(w + b 2 ) 7 2 f3o} + —^-de[sm9( P + b 2 /2)}-cot8^— - =0 (29) 
r sin 6 r 

Eq. (29) shows that typically on a flow expansion time f3e ~ l/7 2 #- Thus for angles 6 > I/7 the lateral 
dynamics is frozen-out for ultra-relativistic flow, and the geometry of the outflow is likely to be determined 
in the acceleration region near the light cylinder. 

One particular stationary outflow configuration contains a quadrupolar current flow concentrated close 
to the polar axis and the equator. This current distribution minimizes the total energy given a total toroidal 
magnetic flux and has been advocated in relativistic stationary winds (Heyvaerts & Norman 2003). In a 
stationary wind, the return current flows at infinity from the pole to the equator; in case of electromagnetic 
explosions, the return current flows along the surface of the bubble - the contact discontinuity. The magnetic 
field in the bubble is inversely proportional to the cylindrical radius. In what follows, we adopt this particular 
current structure. Eventually, the problem of flow acceleration near the source (near the light cylinder of 
the progenitor) and collimation is likely to be solved by numerical simulations. At this moment there are 
no strongly relativistic numerical simulations that trace the flow evolution from the subsonic accelerations 
region, through the special points of the flow to asymptotic infinity. 

There are two important exceptions to the application of these principles. First, as we discussed in 
Section 5.1, at the non-relativistic stage plasma pressure, poloidal magnetic field and tangled component of 
the magnetic field may be important in providing support against hoop stress inside the current regions. 
Secondly, at late times, t > t s , the flow need not be in equilibrium since it changes on dynamical (light 
travel) time scale and, in addition, energy may flow towards the axis where it is dissipated and radiated as 
7-ray emission. 

The minimum size of the core region is the magnetic Debye radius 



^ = fc (30) 

(for a given density n and total current / the relative drift of current-carrying particles in this case is of 
the order of the speed of light). Note that it is orders of magnitude larger than the light cylinder radius 
<~ r s . The corresponding minimal angular size of the core region can be estimated if we use definition (16) 
to eliminate the density 

(31) 

A type of collimation that we envision in case of electromagnetic explosions is somewhat different from 
the conventional "jet" model of AGNs and pulsars. We expect that large Poynting fluxes associated with 
explosive release of ~ 10 51 ergs are sufficient to drive a relativistic outflow over a large solid angle, so that 
during the relativistic stage the resulting cavity is almost spherical. But the Lorentz factor T of the CD may 
be a strong function of the polar angle. 
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5.6. Quasistationary force-free outflow, r < ct s 

At small radii, r p h < r < ct s , the interaction with the circumstellar matter is not important, so that 
the flow may be considered as stationary. At this stage the poloidal magnetic field is small, flow is ultra- 
relativistic and moving ballistically. We can drive general equations governing the behavior of such flows. 
Although the flow may still be supersonic, so that matter inertia is important, in the asymptotic regions 
most important forces normal to the field are electromagnetic (Nitta 1995; Heyvaerts & Norman 2003). 

Consider a simplified model problem of steady state force-free wind carrying only toroidal magnetic 
field. The time-independent force-free flows are described by the stationary Maxwell equations with time 
derivatives set to zero: 



curlE = 

curlB = j (32) 

(factor 47r has been absorbed into definition of current density). From Eq. (32), it follows that the electric 
field is a gradient of a potential, E = — V$, and 

sm9 

= a r $A$ (33) 

r 

Eliminating A<I> from cqns (33) we find 

gggftg- »™gg frg = o (34) 
r smfc> 

Which means that electric potential is a function of total current enclosed within a magnetic loop at polar 
angle 9: 

$ = $(/) (35) 
The equation for the current I(r, 9) then becomes 

+ AI$' Z + (V/) 2 $'$" = (36) 



sin 2 Or 2 

This equation resembles closely the relativistic Grad-Shafranov equation {e.g. Beskin 1997). The main 
difference is that under assumption of zero poloidal field the conventional flux function loses its meaning. 
In the standard theory of the Grad-Shafranov equation, there is a poloidal field which is related to the 
derivative of the flux function. We replace the flux function with the current enclosed by a given field loop 
and treat this as an independent variable. Equation (36) is considerably simpler than the full relativistic 
Grad-Shafranov equation. It has one free function - distribution of electric potential that in the steady state 
determines uniquely the current distribution. In the steady state the plasma experiences radial and lateral 
drift with velocities 

(3= § = -$'sin6»a e lnV7, A = = $V sin 9d r In VI (37) 
B B 

For radial motion A = 0, I = 1(9), equation (36) gives 



de{l(9)^T~P) =0 



(38) 
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which can be integrated 

J(fi)7r^=/o/l^ (39) 

where Jo is the axial current and flo is the ratio of the axial charge density per unit length to axial current. 
Note that condition (3 < 1 requires 7 7^ 0. Thus, subluminal expansion along conical surfaces requires 
presence of a line current (see also Heyvaerts & Norman 2003). Further solutions of this equation can be 
written down but will not be considered here. 



6. Electromagnetic Formalism 

6.1. Relativistic Force-free Electro-dynamics 

The conventional method for handling relativistic, magnetized flows is to use the relativistic extension 
of regular, non-relativistic magnetohydrodynamics. Relativistic MHD (RMHD) is considerably more com- 
plicated. In the non-relativistic case the displacement current is neglected, <9 t E — > 0, and a charge neutrality 
is assumed, p e — > 0. Both of these assumptions may be violated in relativistic plasmas. In addition, when 
calculating dynamic properties one needs to take into account the inertia of the electromagnetic field. How- 
ever, there is a simpler extension, the relativistic force-free approximation (RFF), which is appropriate when 
the plasma is sufficiently tenuous and subsonically moving, so that its inertia can be ignored. On the other 
hand, the local microscopic plasma time scale (e.g. plasma frequency) typically is much shorter than the 
global dynamical time scales and there is plenty of charge available to screen the component of electric field 
along the magnetic field. In this case we can neglect the inertia of the plasma particles but have to include 
their electromagnetic interaction. This is relativistic force-free (RFF) approximation. (Inertial contributions 
may be included as perturbations to magnetic forces.) 

There are two equivalent ways of deriving RFF equations. First, they can be derived from the two 
Maxwell equations (factor An has been absorbed into definitions of currents and charge densities): 

BE 

- VxB-j (40) 

- = -VxE (41) 

with the current density perpendicular to the local magnetic field determined by the force-free condition, 

pE+jxB = (42) 

This immediately implies that the invariant E • B = and its temporal derivative can be set to zero; 
in addition, electromagnetic energy is conserved, E • j = 0. This allows one to relate the current to the 
electro-magnetic fields 

. _ (E x B)V - E+(B-VxB-E-Vx E)B 
J - ^2 

This may be considered as the Ohm's law for relativistic force-free electrodynamics. 



(43) 



For consistency of RFF, one also needs to assume that the second electromagnetic invariant is positive, 
B 2 — E 2 > 0. This implies that there is a reference frame where the electric field vanishes Equivalently, the 
electromagnetic stress-energy tensor can be diagonalized. Note, that there is no mathematical constraint 
that would ensure that B 2 — E 2 remains positive. Thus, even if initially B > E everywhere, both numerical 
and analytical solution (e.g. Appendix F) may have regions where this condition becomes violated. The 
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implies that the RFF is no longer valid, so that either plasma inertia or dissipation should be taken into 
account (Section 6.2). 

The conditions E • B = and B 2 > E 2 allow us to define an electromagnetic velocity v = E x H/B 2 
perpendicular to the magnetic field. The velocity along the field is not defined. This degeneracy comes 
from a neglect of plasma dynamics associated with non-force-free motion of plasma along the field (however 
strong the magnetic field is, it does not influence the plasma motion along the field). Note, that under the 
force-free approximation all plasma species drift across the magnetic field with the same velocity. Thus, the 
current perpendicular to the magnetic field is exclusively due to the plasma charge density. 

RFF equations (40-43) represents a simple evolutionary dynamical system (Uchida 1997a; Komissarov 
2002), which can be solved numerically (e.g. Komissarov 2001). When one includes the constraints E • B = 
V • B = 0, there are four independent electromagnetic variables to evolve and four characteristics along 
which information is propagated. In the linear approximation, these correspond to forward and backward 
propagating fast and intermediate wave modes with phase speeds c and ck • B respectively (k and B are 
corresponding unit vectors). 

RFF dynamics can be developed in a manner that is quite analogous to regular hydrodynamics, with 
the anisotropic Maxwell stress tensor taking the place of the regular pressure and the electromagnetic energy 
density playing the role of inertia (cf. Uchida 1997a; Komissarov 2002). There is an important difference, 
though, in that the existence of a luminal fast mode means that electromagnetic "flows" do not become truly 
"supersonic" (see Section 6.2). 

We would like to stress an important difference between non-relativistic and relativistic force-free fields: 
the relativistic force-free theory is dynamic. In laboratory (non-relativistic) plasma the notion of force-free 
fields is often related to the stationary configuration attained asymptotically by a system (subject to some 
boundary conditions and some constraints, e.g. conservation of hclicity). This equilibrium is attained on 
time scales of the order of the Alfven crossing times. In strongly magnetized relativistic plasma the Alfven 
speed may become of the order of the speed of light c, so that crossing time becomes of the order of the light 
travel time. But if plasma is moving relativistically its state is changing on the same time scale. 



Force-free electrodynamics assumes that the inertia of plasma is negligible. This approximation is bound 
to break down for very large effective plasma velocities, when electric field becomes too close in value to 
magnetic field, E — > B. The condition that inertia is negligible is equivalent to the condition that the effective 
plasma four- velocity, u ~ E/(By/l — (E/B) 2 ), is smaller that the Alfven four- velocity in the plasma. In 
relativistic plasma 



6.2. Applicability of force-free approximation 




(44) 



This puts an upper limit on the value of electric fields consistent with force-free approximation 



B 2 -E 2 
~~E 2 



1 

> - 



a 



(45) 



A second constraint comes from the fact that in RFF the charge number density n e may be comparable 
to the total plasma number density n and the current density may reach j ~ nec. Force-free electrodynamics 
assumes that there are enough charge carriers in the plasma to assure that the condition E B = is satisfied. 
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This condition can be violated if the plasma density is too small and/or if fields change on very small scales 
(e.g. due to development of electromagnetic turbulent cascade which will bring energy to smaller and smaller 
scales). In this case plasma becomes charge-starved. To estimate charge-starved condition, assume that a 
typical fluctuating amplitude of magnetic field at the scale I' is by (both quantities measured in the plasma 
rest frame). The corresponding current j[, <~ cbii/Air is limited by the number density of real charges n': 
j[, < 2en'c, which gives (Thompson & Blaes 1998) 

If this condition is violated, the electromagnetic perturbations behave more like electromagnetic waves than 
MHD waves. In particular, for smaller densities the condition E • B = is not satisfied, so that electric fields 
can accelerate particles, which leads to dissipation of electromagnetic energy. This can be an important 
acceleration mechanism. 

Thirdly, non-ideal effects (such as resistivity) may lead to violation of the ideal RFF approximation. 
It is possible to formulate equations of resistive RFF by introducing a macroscopic resistivity (Lyutikov 
2003a). In resistive RFF, the dynamics is still controlled by the electromagnetic field, but now the currents 
that support these fields may become dissipative. Resistive effects in magnetically-dominated plasma are 
somewhat different from the non-relativistic analog. In a force-free plasma, conduction currents only flow 
along the magnetic field in the plasma rest frame, so that only the component of current along the field is 
subject to resistive dissipation. Ohm's law in resistive RFF becomes (Lyutikov 2003a) 

(j ■ B) + d t ^ v /l-((ExB)/B 2 ) 2 (j • B)) = (B • V x B - E • V x E)B (47) 

In spite of the differences between resistive RFF and non-relativistic plasmas, effects similar to familiar 
resistive instabilities, like tearing mode, should develop in RFF (Lyutikov 2003a). This has important 
implications for the numerical modeling of RFF fields: similar to fluid dynamics, where sound waves turn 
into shocks, RFF fields tend to form small scale structures in which resistive effects becomes important. 



6.3. Relativistic force-free versus MHD 

It is instructive to contrast this electromagnetic approach with the relativistic MHD (RMHD) formalism 
that is currently used by most investigators (e.g. Phinney 1982; Camenzind 1986; Takahashi et al. 1990; 
Takahashi 2000; Park 2002). In the RMHD formulation dynamic equations include incrtial terms and pressure 
gradients. This means that a fluid velocity field must be tracked, along with enthalpy density and pressure. 
The electric field is supposed to be related to current and magnetic fields through an Ohm's law. Under 
relativistic MHD, the constitutive relation, Eq. (43) must be augmented with inertial terms. The equations 
are still evolutionary, though more complex - there are seven independent variables to evolve, with seven 
independent wave modes (two fast, two intermediate, two slow and one adiabatic) to follow. 

The introduction of these extra complications, when the inertia of the plasma is relatively small, can 
be questioned on several grounds. Firstly, it is assumed that the electric field vanishes in the center of 
momentum frame of the plasma. This is not guaranteed by plasma physics in a tenuous, relativistic plasma. 
Secondly, it is usually assumed that the plasma slides without friction along the magnetic field. In other 
words, there is no dissipation. This is unlikely to be the case in the face of instabilities and radiative drag 
(Bcskin & Pariev 1993). Thirdly, it is generally supposed that plasma is conserved. This is untrue in GRBs 
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where pair creation may be important until far out in the outflow. Finally, there is the common assumption 
that the particle pressure tensor is isotropic. In practice this is rarely the case in observed plasmas. In 
contrast to that, the magnetic flux is a much better conserved quantity. 

The relativistic force-free equations may be derived from relativistic MHD formulation in the limit of 
negligible inertia. This offers an advantage that the system of equations may be set in form of conservation 
laws and can be easily generalized to general relativistic form. The basic equations of RMHD include the 
relativistic Ohm's law and the relativistic dynamics, Maxwell's and mass conservation equations Lichnerowicz 
(1967); Uchida (1997a); Komissarov (2001): 

= 0, (48) 
F*f v = 0, (49) 

(pu");m« = (50) 

where the stress energy tensor is a sum of contributions from matter and fields: 

~ (m) (em) 

1 b 2 

= K F Va ~ ^F afj F af) g^ = b 2 u»u v + —g^ - b^ 1 b u (51) 

where w is the plasma proper enthalpy, p is proper plasma density and p is pressure, b a = \r\a.p i iv'uPF ilv 
are the four-vector of magnetic field, Levi-Civita tensor and electro-magnetic field tensor, u a = (7,7/?) 
are the plasma four-velocity, F^ v is electromagnetic field tensor. An additional constraint is that the first 
electromagnetic invariant is zero: 

F^v *F^ = (52) 

where *F fiV is a dual electromagnetic tensor. 

It has been also implicitly assumed in the derivation of these equations that the second electromagnetic 
invariant F^ v F^ v ^ 0, so that the electromagnetic stress energy tensor can be diagonalized. This assumption 
is important since we are interested in the limit when matter contributions to the stress energy tensor vanish; 
the possibility of diagonalization of the electromagnetic stress energy tensor distinguishes MHD and vacuum 
electromagnetic fields where such diagonalization is not possible. 

When we set the matter contributions to the stress energy tensor to zero, w,p — > 0, the RMHD equations 
become equivalent to the RFF equations with the exception that the velocity along the field is not defined 
anymore. This can also be checked by combining RFF equations in the form of conservation laws 

^(£ 2 +B 2 ) + V-(ExB) = 

<9 t ExB = VxBxB + VxExE+(V- E)E. (53) 



7. Relativistic expansion of the magnetic shell, r > ct s 

After breakout, the electromagnetic flow is accelerated by pressure and magnetic forces and becomes 
super-fast-magnetosonic. At the photosphere the lepto-photonic component decouples from the magnetic 
field, so that a increases to a ~ 10 8 — 10 10 , so that the wind becomes sub-fastmagnetosonic again. Beyond 
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this point the outflow is strongly magnetically dominated and can be described by RFF equations. The wind 
is then further accelerated by magnetic pressure, while its expansion is determined by the interaction with 
external medium. In this section we study the dynamics of such magnetically-dominated explosions. 

There are two principal differences between this work and other (mostly stationary MHD-type) ap- 
proaches. First, the flow is commonly assumed to expand freely without interaction with the circumstellar 
medium, so that the flow dynamics is determined by the internal structure of the flow. MHD flows usually 
cross fast magnctosonic critical surface after which they become causally disconnected from their source (Gol- 
drcich & Julian 1970). A condition that the flow remains regular at the critical surface is what determines 
the global properties of the wind. Unlike MHD, force-free flows are sub-fastmagnetosonic so no conditions 
at the fast critical surface appear. In this case it is the interaction with boundaries that determines the 
properties of the flow (similar to subsonic hydrodynamic flows). In the case of freely expanding magnetized 
wind the terminal Lorentz factor Too depends very sensitively on details of geometry: for radially expanding 
MHD flows Too ~ Michel (1969); Goldreich & Julian (1970). Thus, the distinctive feature between MHD 
and force-free flows is whether the wind becomes fast supersonic (MHD regime, T > ^fa) or not (force-free 
regime, V < \fa). In case of GRB, beyond the photosphere a <~ 10 9 , so it the interaction of the force-free 
wind with the external medium (and not the internal wind dynamics) that is likely to limit acceleration to 
r - 10 4 (see below). 

Secondly, the flows under consideration are strongly non-stationary, so that their dynamics is qualita- 
tively different from stationary solutions and is more reminiscent of the strong explosion problem (Sedov 
1969). Relativistic effects introduce an important difference, though. Primarily, the differences lie in the 
causal structure of the shell. It takes a very long time for a wave to be reflected by the relativistically 
expanding surface of the shell (contact discontinuity) and return to the origin. This is generally true of 
ultra-relativistic flows so that stationary solutions, which take a long time to be established, can be quite 
misleading. 

As we have discussed above the long time scales for establishing a "feed-back" on the source have an 
important consequence for energy dissipation in the flow: there is no necessity to destroy magnetic flux 
through ohmic dissipation, until the wave can actually propagate back to the source. Stated another way, 
there need be little resistance in the electrical circuit. The effective load can consist of the performance of 
work on the expanding blast wave. This is where most of the power that is generated by the central magnetic 
rotator ends up. Thus, until the outflow becomes non-relativistic, the distinction between the inertial load 
and a different, dissipative load is quite unimportant for the behavior of the central source. As long as 
the expansion remains ultra-relativistic, it is a very good approximation to impose a Sommerfeld radiation 
condition on the solution at the source and to ignore the reflected wave. 

The general problem of relativistic force-free expansion is quite complicated: we need to solve for 
the time-dependent, three-dimensional evolution of electric and magnetic fields with a non- linear coupling 
through the force- free Ohm's law (43). (note that even normal modes in such plasma experience strong 
refractive phenomena). To make the problem tractable we make several simplifying assumptions. First, we 
limit our approach to distances much larger than the source scale r s and, since at large distances magnetic 
field is dominated by the toroidal fields, we neglect poloidal magnetic fields altogether (thus we also neglect 
the angular momentum loss by the source). Secondly, instead of treating the flow acceleration and collimation 
at the source, we assume that the central source can be represented as an inner boundary condition through 
which a toroidal magnetic flux is injected with a given total luminosity and a given polar angle dependence 
L{6). 
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We assume also that at large distances, r > r p h, the outflow may be described by ideal RFF equations. 
In this case the expansion of the relativistic shell will be subsonic, and virtually ballistic (radial). The energy 
input of the central source then can be represented as fast modes propagating from the source to the contact 
discontinuity. Motion of the CD is determined by the pressure balance between the Poynting flux from the 
source and the ram pressure of the ISM and thus depends on the source luminosity La(t') per unit solid 
angle <Kl at the retarded time t' such that R(t) = t — t' . Typically R(t)/c ~ r 2 t sourco . 

A particular form of the forward-propagating wave train (its temporal variations, amplitude and lateral 
distribution of the energy flux) depends on the history of the source activity and detailed properties of the 
source. Since the central source is expected to vary considerably during its activity, the pressure on the 
inner boundary of the CD will be fluctuating as the new pulses from the source start to contribute to the 
pressure. This will be reflected in the "jitter" of the CD motion. At later times multiple reflections from 
the contact discontinuity and the center become important as well. As we have discussed in Section 5.5 the 
lateral distribution of the magnetic field and energy fluxes in the outgoing waves is likely to consist of a line 
current. Then, a very important observational consequence follows. For outgoing waves corresponding to 
the line current the distribution of Poynting flux is L ~ 1/sin 2 9, implying that the central source releases 
an equal amount of energy per decade of 6. This would later translate into a similar form of the forward 
shock. 

The form of the ingoing reflected wave depends both on the outgoing wave and the form of the contact 
discontinuity. Since in the strongly relativistic regime the wave will be reflecting from receding CD the 
amplitude of the reflected wave will be ~ T 2 times smaller. 

The relativistic stage of the shell may be separated into two stages: "coasting" and "self-similar" 
depending on whether or not most of the fast waves emitted by the central source have caught up with the 
CD. 

• "Coasting" electromagnetic shell (ct s < r < r s h = (Lat^/pc) 1 / 4 ) (in the observer frame this phase 
lasts also ~ 100 s). At r > r p h the shell becomes a relativistically expanding shell of thickness 
<~ ct s <~ 3 x 10 12 t s ,2 cm. The shell still contains toroidal magnetic field but the current now detaches 
from the source and completes along the shell's inner surface. At this stage the CD is constantly re- 
energized by the fast-magnetosonic waves propagating from the central source. The average motion of 
the CD is determined by the average luminosity, V ~ (Ln/pc 3 ) 1 / 4 /- -1 / 2 (in a constant density medium) 
or T ~ const (in a p ~ r~ 2 wind). In addition, there is a "jitter" of the CD, in response to the "jitter" 
in the source luminosity. The internal structure of the magnetic shell is a messy mixture of the outgoing 
waves from the source and the ingoing waves reflected from the CD, similar to a pre-Sedov phase in 
hydrodynamical explosions. Unlike the case of a hydrodynamic blast wave with energy supply, no 
internal discontinuities form inside magnetic shell. 

• Self-similar electromagnetic shell (r s h < r < tnr = (Lat s / pc 2 ) 1 ^). After one dynamical time scale, 
td ~ T 2 t s ~ 3 x 10 16 cm, all the regions of the shell come into causal contact - most of the waves 
reflected from the CD have propagated throughout the shell. As the expanding shell performs work 
on the surrounding medium its total energy decreases; the amount of energy that remains in the shell 
during the self-similar stage is small, ~ i?o/r 2 . Most of the energy is still concentrated in a thin shell 
with AR <~ R/T 2 near the surface of the shell which is moving according to T <~ \J E^j pc 2 r^l 2 (in a 
constant density medium), or T ~ r^ 1 / 2 (h a p ~ r~ 2 wind). The energy of the shell decreases only 
weakly with radius, dE/dt ~ 1/r in constant density and dE/dt ~ 1/r 3 in a wind, so that the surface 
of the shell keeps moving relativistically as long as the preceding shock wave is moving relativistically, 
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until r ~ (E n /pc 2 ) 1/3 ~ 10 18 cm - the shock never becomes completely free of the shell. Interestingly, 
the structure of the magnetic shell (in particular the distribution of energy) resembles at this stage 
the structure of the hydrodynamical relativistic blast wave (Blandford & McKee 1976) (Section 8.1.1, 
m=3). 



7.1. RFF equations for axisymmetric asymptotic flows 

In the asymptotic domain r>r s the ejecta will be dominated by toroidal magnetic field and, since we 
assumed that all small scales field variations are smoothed out, it will be axially symmetric. In this case the 
only nonzero components of the fields are B = B^, E ri Eg and the RFF equations (40) and (43) give 

8 t B = --d r {rE g ) + -d E r 
r r 

d t E = --d r (rB) + ^ ( —L-dg(sm 8E e ) + ^d r (r 2 E r ) 
r B \rsm0 r z 

1 Q Eg ( 1 Q , . „,. , , I ., , , 



d t E r = -_-d e (sin0B) - ——de(sm0E e ) + -^d r {r z E r ) (54) 

The charge density associated with a given solution is given by 

p=^d r {r 2 E r ) + -^J-E S m6 (55) 
r z r sin o0 



and the currents are 



E e _ Eg ( 1 2 , . 1 8 



B B \r z rsmvdv 



30 = -p-£ = ~ -d r (r 2 E r ) + ———Esm0 (56) 
B B \r z r sin o 06 



After some rearrangement the system (54) gives 



-d t (E 2 + E 2 + B 2 ) + -^d r r 2 E B B — dg sm9E r B = 

2 ' r z rsm.0 

d t E B + \d r r 2 (B 2 + Ef) ]— dg sin6E r E e - \d r r 2 E r = 



E E E B E 

d t E r B + —d r rE + -4rd T r 2 E r + e —d e sm6E e -dgsuvOB -d e E r = (57) 

r H rsmV rsmo r 

To connect to RMHD description we note that the plasma drift velocity is (Eg/B)e r — (E r /B)eg. Setting 
B = ■yb, Eg = (iyb, E r = X-yb and the system (57) gives eqns (C6). 

Equations (54) are the main equations that determine dynamics of the expanding magnetic shell. 



7.2. Boundary condition on the contact discontinuity 

To find the internal structure of the expanding magnetized shell the dynamical equations should be 
supplemented by boundary conditions on the contact discontinuity separating the relativistically expanding 
magnetic shell and external medium. The first boundary condition requires that the pressure of the magnetic 
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field of the ejecta should be balanced by the pressure of the external medium. The pressure of the external 
medium is the ram pressure if the forward shock has not formed, or the kinetic pressure of the shocked 
material if the forward shock has formed. The two cases differ only by a numerical coefficient n (appendix 
B). 

As the shell propagates into the surrounding medium it reflects particles with typical Lorentz factor 
<~ 2L 2 . In the frame of the shell there must be a pressure equilibrium on the contact discontinuity. The 
magnetic pressure inside the shell, -B 2 /(87iT 2 ) should balance the ram pressure of the incoming external 
medium , 2KL 2 /9 cxt c 2 , (here p ext is the density of the external medium and k is a coefficient of the order of 
unity relating the ram pressure before the shock and kinetic pressure on the CD): 

2 K r 2 /3 2 Poxt c 2 = ^ (58) 

where B is the shell magnetic field measured in the lab frame. Thus, the magnetic field in the plasma rest 
frame should scale as b ~ T, while magnetic field in the lab frame scales as B ~ Tb ~ T 2 . 

Alternatively, a force balance on the CD can be considered as a Lorentz force of magnetic field and 
the surface current acting against the ram pressure. Consider, for example, magnetic field occupying a half 
space x > and moving with velocity (3 x e x + /3 z e z , so that the fields are given by B = B Q(x — (3 x t)e y , E = 
—(3 x BqQ{x — (it)e z + [3 z BqQ(x — (3t)e x where 8 is the Heaviside function. Then from Maxwell equations 
the surface current g is 

g=^5(x-(3t)e z , (59) 
where L 2 1/(1 - (/3 2 + /3 2 )), so that the Lorentz force on the surface is g x B = Bq/AttT 2 . 

The second boundary condition comes from the requirement that the flow should smoothly connect to 
the CD, which implies that in the frame of the CD Eg on the CD equals zero, or, alternatively, the component 
of the electromagnetic velocity normal to the surface of the CD should equal the velocity of the CD. If at a 
given point the CD has a velocity v CD = {cos a C D, sina C£ >, 0} in spherical coordinates (r, #,</>), then 

E e cos a C D - E r sin a C D 

2 = v C D (60) 



where tanacc = dg\nR(9,t) <~ 0(1/P 2 ) and vqd = \/l — 1/r 2 . Note, that in order to calculate the angle 
olcd between the radial direction and normal to the CD one needs to know the whole history of the shell 
expansion since it depends on the position of the CD at a given time R(9, t). On the other hand the velocity 
of the CD vcd depends only on the source luminosity at a retarded time. 

The third boundary condition, for the radial electric field E r (R), determines the surface charge density 
on the shell. 



7.3. Motion of the contact discontinuity 

Next wc determine an average motion of the CD given an average source luminosity. At the coasting 
phase, equating the source luminosity Ln(t') with the Poynting flux at the CD Ln(t') = B 2 R{t) 2 (3 / Air and 
using boundary condition (58) we find 



L n (f) = ^pc 3 T 4 R(t) 2 f3 3 



(61) 



-30- 



If during its activity period the source produces an approximately constant luminosity Lq_ , then for a constant 
density medium p(r) = p 

T CD = (io/4 K p c 3 ) 1 / 4 i?- 1 / 2 , (62) 
while for a wind environment p(r) = po(ro/r) 2 

T CD = (L n /A K p r 2 c 3 ) 1 / i = const (63) 

If the distribution of luminosity is as argued above, Lq = L a 9~ 2 , then Tcd oc l/\/smd (see Section 8.2). 

The work that the shell makes on the circumstellar medium per unit time and unit solid angle is (see 
Eq. (59)) 

pdV/dt = g x B\ CD R 2 (3 = Aupc^T 2 R(t) 2 [3 (64) 

Thus, in case of constant density the most efficient energy transfer to the forward shock occurs at the end 
of the coasting phase, pdV/dt oc R, R < 2T 2 t s c, while for the wind environment pdV/dt <~ const. 

By the end of coasting phase most of the waves emitted by the central source will have reflected from the 
CD, reached the center and reflected again. At this stage the structure of the shell will become self-similar. 
At the self-similar stage T <~ i?~ 3 / 2 for constant density and T <~ i? -1 / 2 for the wind environment. By the 
end of the coasting phase, a large fraction of the energy of the shell will have been transferred to the forward 
shock. The transfer of the remaining small fraction (of the order of E tot /T c , where r c ~ 100 is the Lorentz 
factor at the end of the coasting stage) to the forward shock becomes very inefficient (pdV/ dtdQ <~ 1/r for 
both p = const and p oc 1/r 2 ), so that the energy inside the shell decreases only logarithmically with time 
until the motion becomes non-rclativistic. 



8. Internal structure of the magnetic shell 

To determine the internal structure of the magnetic shell one needs to solve Eqns. (54) together with 
boundary conditions (58-60). This is a quite complicated system and should generally be solved numerically 
In this section we first consider an idealized case of purely radial motion which allow exact solutions, and 
then derive equations and obtain approximate solutions for our preferred case, when the outgoing wave 
consists of a line current. 



8.1. Radially expanding magnetic shell 

To illustrate the behavior of the solutions we consider in details a particular case of radial propagation 

of both outgoing and reflected waves. In this case exact solutions can be found. For radial propagation 
E r = and the relevant equations become 

■at = -d~r rE (65) 

^ = _IA rB ( 66 ) 
dt r dr 

dg (sin 2 6(B 2 - E 2 )) = (67) 



Eqns (65-66) can be reduced to a 1-D wave equation, which solutions can be expressed as a superposition 
of simple wave (Landau & Lifshits 1975a). Taking into account that in the force-free plasma the fast 
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magnetosonic waves propagating across magnetic field have velocity c, the simple waves are given by r = 
±t + f(v), where f(v) is some function determined by the boundary or initial conditions. Allowing for 
angle-dependent velocity of expansion the general solution of (67) becomes 12 



B = 

E = 
(B 2 



h{t-r)z{9) + 



fi(t-r)z(0) 



Mt + r) 

z{9) 
h(t + r) 



z{6) 



9(0) 
9(0) 



£ 2 )sin 2 



h{r,t) 



(70) 

(71) 
(72) 



The two terms in each expression are fast modes propagating outward and inward. The form of the functions 
fi and /2 is determined by the boundary conditions at the surface of the shell. 13 



For given fields (72) the radial velocity of the plasma is defined as 

_ hzjef - f 2 

' hz(6? + h 

while the charge density p — divE and current j = curlB — 9 t E become 



(73) 



P 

jr = 



d e (sm0z(0)g(0)) . d e (sm0g(6)/z(9)) 



sin Br 1 
dg(sm9z(8)g(8)) 

sin Or 2 



h 

fi + 



sin 9r 2 
d e (sm9g(9)/z(9)) 



sm 



9r 2 



h 
h 



(74) 



In the plasma rest frame electric field is while magnetic field is b — 2\J fifig(0)/r. (Note that it depends 
both on the amplitude of the forward and backward waves.) 



The boundary conditions (58-60) become 



(F lZ (8)+F 2 /z(8)f 

r 2 



(F lZ {0) + F 2 /z(9)j 2 = alt 2 sin 2 9T 4 f3 2 
and w = 167TK/9 ext c 2 and F = f(r = R(t)). We can solve the system (75-76): 



Fi = 



l-P2z{9)g{9) z{9)g{8) 



(75) 
(76) 

(77) 



12 It is instructive to derive these equations using the RMHD formulation. Under assumption of radial motion Eqns (C6) 
become 

/3~/d t b + bd r ~/ = 

^drW + = Q (6g) 

r 

Which using substitution rb —* cxp{6i}, 7 — > coshj;, /3 — ► tanhj; reduces to 

d t bi + d r y = 

9rbi + dty = (69) 
with a general solution 61 = c 1 (6)(f 1 (t - r) + h(t + r)) + c 2 (0), y = ci(6»)(/i(t - r) + f 2 (t + r)) + c 3 (6»), which reduces to (72). 

13 In case of a dominant toroidal magnetic field, these equations can be generalized for MHD plasma, when inertial and 
pressure effects are important. For propagation perpendicular to magnetic field lines (and only in this case) simple waves in 
MHD may be reduced to hydrodynamic simple wave with a special form for the equation of state, including magnetic field. 
Then the one-dimensional rclativistic simple waves become r = t(v ± Uf)/(1 — vuf) + f(v) where uj is the fast magnetosonic 
speed (Landau & Lifshits (1975a)). 
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_ R^z(9) _ R^z{9) 
2 1 + /3 2. 9 (0) ~ Ag{9) 

where approximations assumes r> 1. In addition, Eq. (72) gives 

9(0) = ~?-z (79) 
sin 9 

Equations (78-79) define possible angular dependence of the solution since both F\ and F2 are angle inde- 
pendent by definition. 

Since for strongly relativistic motion R ~ t. from Eqns. (77-78) it follow that there are two possible 
choices. First, w = const = w , z{9) = g{9) = 1/ sin0, T oc g(9) = T (t)/ sm6 with the fields given by 

B _ fi(t-r) | h{t + r) 
r sin 2 6 r 

E= h^_MttA (80) 

rsin 2 6» r 

The second possible choice is W oc g{9) 2 =wq/ sin 2 9, z(9) = const, T(0) = const = T (t). 

h(t-r)+f 2 {t + r) 



B = 



rsmf 



e, . *('->:*('+-> (81) 

r sine/ 



In both cases the functions /1 are and fi are found from boundary conditions 



F 1 = t^T (ty (82) 



F2 = (83) 

Note, that the amplitude of the reflected wave /2 is Tq times smaller than that of the forward propagating 
waves fi due to the Doppler shift during reflection off the relativistically moving surface of the shell. Eqs. 
(82-83) describe the jitter of the CD as a function of the source history of activity through the dependence 
of fi on the retarded time (this is true only at the coasting stage, when no waves reflected from the CD have 
reflected from the origin and reached the CD for the second time). 

In the ultra-relativistic case r > 1 a further simplifications is possible. Since for relativistic motion 
R(t) ~ t the function / 2 can be found explicitly 



while /1 can be found for any given expansion velocity of the shell T (t) and R(t). 



8.1.1. Self-similar solutions 

A particularly simple and analytically tractable case of the radial expansion of the magnetic shell is 
when the Lorentz factor of the contact discontinuity is a power- law function of time T 2 = r 2 n (£/£ in )~ m . 



-33 - 



This includes a point-like explosion and a power-law luminosity variation of the central source. We will call 
these self-similar solutions. In this case R(t) = t (l — l/(2(m+ 1)T 2 )), so that at the boundary 

Fi ~Fi(t/(2(m + l)r 2 )) 

F 2 ~F 2 (2t) (85) 

We can then find functions /i and f 2 

h{z) = v^r 2 „ (2(m + i)r 2 „z) (1 - m)/(1+m) = 0^ 2 * 2ro/(1+m) (2(m + i)r 2 z) (1 - m)/(1+m) 

/ 2 (*) = \A^/8 (86) 
Equations (86) determine the amplitudes of the outgoing and the ingoing waves for self-similar solution. 

We can find the structure of fields in the shell for the self-similar expansion. Keeping only the leading 
terms in F 2 we find 

B oc fl{t ~ T) = V^r 2 - [2(m + l)r 2 (l - r/t)] = ^ 2 * { i- m)/( i +m ) (g7) 

where x = 2(m + l)r 2 (l - r/t) (c.f. Blandford & McKee 1976) (see Fig. 6). For a given time t the 
magnetic field depends on a slowly varying function 1/r multiplied by fast varying (on a scale i?/r 2 ) variable 
^(i-m)/(i+m)_ p or exam pi e; f or a constant energy source in constant density medium, m = 1, B oc y/woT 2 n /r 
(times some angular dependence), which is constant in time inside a shell. 

Result (87) implies that in the self-similar case the energy density of the magnetic field near the surface of 
the CD has the same distribution as the kinetic energy density in the case of relativistic shock wave Blandford 
& McKee (1976). This is a surprising result, since the two systems are completely different (electromagnetic 
shell and hydrodynamical shock wave) and solutions come from different equations (Maxwell and Euler). 
Generally, for 1 < m < 3 the magnetic field in the self-similar solutions piles up near the contact discontinuity, 
concentrated in a layer of thickness R/T 2 . In this case the above solutions may be obtained by the Blandford- 
McKee method (which assumes that the most of the energy is concentrated near the outer boundary of the 
cavity) from the RMHD equations (see appendix C) . To continue our analogy with the hydrodynamic case, 
we note that the above force-free solutions may be interpreted in terms of RMHD if we define a Lorentz 
factor 

2 = & = ll+fj + Vih _ Jj_ = t 2 2 (1 _ m)/(1+m) , . 

7 5 2 -£ 2 4/i/a 4/ 2 r (r + t) X [ ' 

Note that near the CD, i~rwe get 

llearCD=T 2 X {1 - m)/{1+m \ (89) 

again similar to the case of hydrodynamical blast wave (Blandford & McKee 1976). 

Two kinds of self-similar solutions are most interesting. The first is the m = 1 case corresponding 
to a constant luminosity source and expansion into a constant density medium. In this case fi(t — r) = 
woT 2 n ti n — const and the fields are given by 



r V^o ( r 2 t + t + r \ w V^o ( r 2 t t + r 

B= V^(rlM t + r\ V^(r]^^t + r\ 

r \ sin 2 6 8 / r \ sin 2 6 8 J v ' 
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The second is the m — 3 case corresponds to the late self-similar stage when most of the energy of 
magnetic shell has been transferred to the forward shock. This case becomes applicable when all the waves 
emitted by the central source have been reflected from the CD and reached again the center establishing a 
causal contact. In this case / = V w ol (t 3 ^ 2 ^/x) 7 an d the fields are given by 



Generally self-similar solutions have a very limited application since at the coasting stage they are 
limited by finite activity period of the source, limiting their application to a narrow layer near the CD, while 
at the self-similar stage (when the whole shell comes into a causal contact) the total energy remaining in the 
shell is small - most of the energy has been transfcrcd to the forward shock. 

It is instructive to relate this simple solution to the properties of characteristics. In force-free electro- 
dynamics, there are just two characteristics, a fast mode and an intermediate ( Alfvcn ) mode. The former 
is simply an electromagnetic wave with electric vector parallel to k x B. It propagates with phase and group 
velocity c in all frames (Uchida 1997b; Komissarov 2002). The intermediate mode is best described in the 
shell frame, where the background electric vector vanishes. It has its electrical perturbation perpendicular 
to B' and magnetic perturbation parallel to k! x B' . The phase velocity is k' ■ B' and the group velocity is 
equal in magnitude to c and directed along B' . When we transform to a general frame, the group velocity, 
V g , which is what is important for the transmittal of information, has a magnitude of c and a component 
E x B I B 2 resolved perpendicular to B, with the remaining component, of magnitude (1 — £ ;2 /_B 2 ) 1 / 2 c, 
directed along ±B. The magnetic perturbation is along k x V g . 

The backward-propagating, reflected wave is a fast mode but has a small amplitude when T ^> 1. (In this 
approximation, the intermediate model's group velocity is purely toroidal and carries no radial information.) 
However, the fast mode is able to propagate radially inward and, in effect, create an electromagnetic pressure 
wave which decelerates the electromagnetic velocity of the outflow and reduces it to match that of the contact 
discontinuity. The magnetic field is toroidal and carries information about the current flow. Put another 
way, if we were to change the properties of the load, e.g., by encountering a sudden increase in the ambient 
gas density, which would cause sudden jumps in T, then the boundary conditions on the interior flow would 
change, along with changes in the amplitude of the reflected wave. This allows the interior solution to adjust. 

In this section we have considered, mostly for didactic purposes, two exact solutions for the radial 
expansion of the shell. Both of them arc not physically realizable: solution (80) requires a large (formally 
divergent) axial current and a similarly large distributed return current, while the solution (81) requires an 
unreasonable density distribution. Despite of these limitations these solutions illustrate the principal point 
of the shell dynamics: that the shell motion is a result of the electromagnetic waves reflecting from it, that in 
some cases it is possible to define electromagnetic velocity of the flow so that electromagnetic waves behave 
in a manner similar to hydrodynamics, and that the self-similar structures of an electromagnetic bubble may 
resemble self-similar structure of a fluid shock. 




(91) 
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8.2. Structure of the shell for axial outgoing current 

Next we consider the most relevant case, when the outgoing wave consists of an axial current (and axial 
charge) only. Since generally the reflected wave is non-radial, this case is considerably more complicated, 
so we are forced to make several simplifying assumptions. First, we assume that the outgoing wave carries 
only a line current and no distribute current or charge density. The reasons behind this assumption are 
discussed in Section 5.5. In this case the outgoing wave is radially propagating. Secondly, as we discussed in 
the previous section, the reflected waves are typically T 2 weaker. Thus, we can expand the electromagnetic 
fields in the shell assuming that the reflected wave is weak. Expanding the field in terms of small amplitude 
of the reflected wave e <~ 1/r 2 <C 1, 

B = -%/(t - r) + eB^ 

Ee = ^-J{t - r) + eEf 
r sin 

E r = eE e , (92) 



we find from eqns. (54) 



d t {B^r)+d r {Ef\)-d e E r = Q (!)3) 



d r {B^r) + dt{E { g 1) r) =0 (94) 



sin 

These equations simplify if we introduce 



a t ( s ing(4"- gW)) | Q r + 3,)(E r r») _ o 



Then Eqns. (93-94) may be combined to give 

d ^)-dUe i e 1) )-d r d e e -^=0, (97) 

which shows that generally for E r ^ the inward propagating wave is not radial. 
We can simplify the system (93-95) if we introduce y = (e^ 1 ' — M 1 '): 

r 2 (d t - d r ) y + d e e r = 

d -^ + (d t + d r )y = (98) 
which can be further reduced to a single PDE for one variable y: 

(ft + d r ) r 2 (d t -d r )y + d e (^M) = (99) 

\ sm J 

This non-linear PDE equation describes the dynamics of the reflected wave. 

Equation (99) allows a separation of angular variable. Assuming that y = f(r, t)g(6), equation for g{9) 
becomes 

de (^m +l{l+1)g=0 doc) 

V sm I 
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which has solutions / = Pi (cos 9) where are generalized Legendre polynomials (for I = the solutions is 
/ = d I sin 9 + C 2 cot 9). 

Equation for f(r,t) may, in turn, be simplified if we introduce retarded and advanced time variables 
£i = t + r and £2 = t — r: 

% ((^-6) 2 %/)-«a+i)/ = o (101) 

If a solution /(£i,£ 2 ) to the Eq. (101) is found, then the corresponding fields are found from relations 
e^ 1 ) oc P*(cos#), e r oc de (sin 9P^ ) / sin 9 and 

2%6« = -(%+%)/ 
2%e« = (%-%)/ 

2%e r = -/ (102) 

Solutions of the Eq. (101) cannot be found analytically, except for I = 0, in which case the general 
solution can be written as a sum of two solutions yi = o = yi=o-i + yi=o-2- 

y*=0;i = ^ (ffi(a) + J (/^l^ ) . = (C + In v^ol^) ff 2 (£ 2 ) 

t <-i(; ( «-, lKl)+ ^4-/^) 

*'-K* (6)+ * (ft > + f^ + /<OT*) (103) 

and 

y ; =0;2 = C 2 COt0 ^(a) + J (ftl^ ^) 

er* = \J H^db, bP = -^cot0, e^^^lcote (104) 

where C, Ci and C 2 are constants and if 2 (£ 2 ) and G 2 (£ 2 ) are arbitrary functions (G 2 (£ 2 ) is just a 

correction to the amplitude of the outgoing wave and can be set to zero) . 



Solutions of the Eq. (99) should be matched to boundary conditions. Eq. (58) implies immediately 

T-i r (t) 



Vsin# 



(105) 



In order to use condition (60) we need to know how the angle (Xcd between the radial direction and the 
normal to the CD depends on and how it evolves with time. This can be done only in the self-similar 
regime, when the motion of the CD is a simple given function T (t) — Ti n (t/t in )~ m where T in is the initial 
Lorentz factor at time t in . In this case we find 

cos# . . 

t*na CD = -— ^ (i _ gin0/(2(m - 1)r2)) (106) 

Thus, 

cos 9 , . 

acD » 2^TWo (107) 

which is only weakly dependent on 9 for 9 <C 1. 
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The angle acD ~ 1 /r 2 <C 1 is the angle at which the waves emitted by the central source fall onto the 
CD (acD is measured in the laboratory frame). Note, that the reflected angle, on the other hand, is of the 
order of unity <~ cccd^ 2 . 

Expanding the boundary conditions for o>cd ~ 0(1 /T 2 ) we find 

r (P W - P (1) + E r a CD ) T 2 sin 9 = J Fi (t - R) (108) 

Recall next, that pW — PW cx P^cosfl) and E r cx de(sm9Pl)/ sin#. This implies that in order to 
satisfy Eq (108) the reflected wave should consist of a large number of spherical harmonics P; 1 with temporal 
and radial dependence given by solutions of Eq. (101). 

At this point we stop our analytical approach since the resulting equations are clearly more easily solved 
by numerical simulations. In this section we showed how the dynamics of time-dependent force- free outflows 
should be calculated, gave simple analytically tractable examples and finally wrote down equations describing 
a particular type of electromagnetic explosion, when the outgoing current is confined to a small region around 
the axis. Though we did not solve for the structure of the reflected wave, its properties have little effect 
on the form of the CD since typically the amplitude of the reflected wave is T 2 smaller than of the forward 
wave. For example, the basic property of the expanding shell, dependence of Lorentz factor on the angle Eq 
(105), can be obtained without solving for internal structure. 

The upshot of this section is that our preferred current distribution - with current concentrated near 
the pole, along the CD and in the equatorial region - will produce a nonsphcrical relativistically expanding 
shell. The internal structure of the shell is determined by a combination of outward and inward propagating 
waves and is in many ways similar to hydrodynamic subsonic flow. Our solutions are formally divergent 
on the axis, Yqd oc 1/Vsin#, Lq = L a / sin 2 9 but the total energy flux is only weakly (logarithmically) 
divergent L to t = J L^dVl k, 2ttLq\yiso\9. The outflow then should have a current carrying core. The core 
may, in fact, be cylindrically collimatcd (e.g. Heyvaerts & Norman 2003), but if we assume that its typical 
polar angle is 9o ~ f0 -3 <C 1 (see (31)), then, qualitatively, 

L ^ ~ . L tot « 27TZ-0 In l/0o (109) 

The energy released per solid angle and the total energy of the explosion will have a similar dependence 

S °"^+0 2 ' E tot ^2irE lnl/9 (110) 
This will play an important role in the interpretation of afterglows (see Section 10). 

9. 7-ray emission 

9.1. Dissipation of magnetic energy 

One of the principal implications of the electromagnetic hypothesis is that the conventional model of 
particle acceleration - acceleration at internal shocks - cannot work in this model. In the RFF limit, the fast 
speed is the speed of light so that fast shocks do not form. If we add a limited quantity of plasma and use 
RMHD, the shocks are weak and not likely to be efficient particle accelerators. There are no intermediate 
shocks in the RFF limit, though rotational discontinuities can be present. 
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We therefore propose that the 7-ray-emitting electrons are accelerated by current instabilities during the 
magnetic shell phase. Current-driven instabilities play a major role in a variety of laboratory experiments, 
e.g. TOKAMAK discharges like sawtooth oscillations and major disruptions, (e.g. Kadomtsev 1975), Z- 
pinch collapse, (e.g. Rudakov & Sudan 1997), in the Earth's magneto-tail (e. g. Galeev et al. 1978) and on 
the Sun (e.g. Aschwanden et al. 2002). The development of current instabilities usually results in enhanced 
or anomalous plasma resistivity which can lead to an efficient dissipation of the magnetic field. The magnetic 
energy is converted into heat, plasma bulk motion and, most importantly, into high energy particles, which, 
in turn, are responsible for the production of the prompt 7-ray emission. The conversion of magnetic energy 
into particles may be very efficient. Recent RHESSI observations of the Sun indicate that, in reconnection 
regions, most magnetic energy goes into non-thermal electrons (Bcnz et al. 2003). We propose, that in case 
of GRBs it is the dissipation of magnetic energy that is responsible for particle acceleration. 

Why is magnetic field dissipation is negligible close to the source and become important at large dis- 
tances? We argue that this is because the particle acceleration is suppressed near the central engine by 
efficient pair production which screens out the electric field that led to the particle acceleration. Eventually, 
the optical depth to pair production becomes small so that large electric fields may develop (see Eq. 27). 
This estimate is consistent with the assumption that particle acceleration and emission generation takes 
place at r > 10 16 cm. (If particle acceleration takes place at r 7 _ 7 > 1 all the IC photons that will neces- 
sarily accompany such acceleration will make pairs. This will increase the density of pairs and will shut off 
acceleration. In this regime most dissipated energy goes to pair production, not synchrotron emission.) 

Magnetic dissipation is complicated. After 40 years of intensive research it is not clear at the moment 
what is the correct macroscopic model of magnetic reconnection, Sweet-Parker or Petschek (Biskamp 2000; 
Priest & Forbes 2000). The dissipation of magnetic energy in relativistic, magnetically-dominated plasma 
is likely to be complicated as well. Only the first steps have been made in that direction (e.g. Lyutikov 
& Uzdensky 2003; Hoshino 2002; Larrabee et al. 2002; Lyutikov 2003a). A fundamental problem is that 
in case of magnetic reconnection most of the observed properties depend sensitively on the kinetic details 
of the model. This can be contrasted with the shock acceleration model, where a basic kinetic property - 
spectrum of accelerated particles - can derive using simple macroscopic quantities - shock jump conditions 
(e.g. Blandford & Eichler 1987). 

Overall, the observed radiative properties of magnetic dissipation and internal shocks should be simi- 
lar. Emission generation (synchrotron) is likely to be the same in reconnection models, so that most of the 
well detailed radiation models will still hold (modulo, perhaps, the assumption of the cquipartition mag- 
netic field). On the other hand, the acceleration mechanism is very different, but, unfortunately, hardly 
distinguishable observationally - both internal shocks and reconnection regions represent transient internal 
dissipative regions, which heat the plasma and accelerate particles. This problem is also complicated by the 
fact that, in the emission region, the plasma should be near equipartition (since it emits most efficiently 
then). Unlike with the fireball model, in the electromagnetic model, equipartition is reached by dissipation 
of the magnetically dominant plasma. 

Particle acceleration by dissipative magnetic fields may proceed in a number of ways. The best studied 
non-relativistic example is particle acceleration in reconnection regions either by inductive electric fields, 
resistive electric fields inside the current sheets (e.g. Craig & Litvinenko 2002) or formation of shocks in the 
downstream of reconnection regions (e.g. Blackman & Field 1994). Investigation of the particle acceleration 
in the relativistic regime in only beginning (e.g. Larrabee et al. 2002). Relativistic reconnection may 
also produce power-law spectra of accelerated particles (Hoshino 2002). For example, in the relativistic 
Sweet-Parker reconnection model (Lyutikov & Uzdensky 2003) , if one balances linear acceleration inside the 
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reconnection layer by the resistive electric field, dt£ ~ eEc with the rate of particle escape (proportional to 
relativistic gyro- frequency) , d t \nN(£) ~ w_b(toc 2 /£), one finds 

N(£)^£- 0i " (111) 

where £ is the energy of a particle, N(£) is the particle number and ~ E/B is the inflow velocity. For 
relativistic reconnection the inflow velocity can be relativistic (Lyutikov & Uzdensky 2003). This simple 
estimate is confirmed by numerical simulations (Larrabee et al. 2002). The fact that reconnection models 
can produce spectra which are prohibitively hard for shock acceleration may serve as a distinctive property 
of electromagnetic models. 

In addition to acceleration at reconnection layers particle acceleration may occur through formation 
of a spectral cascade of nonlinear waves in force-free plasma which transfer energy to progressively larger 
wave vectors until this energy is taken up in accelerating a population of relativistic electrons and positrons. 
Formation of turbulent cascade may be initiated, for example, due to development of dynamic instabilities 
of the CD (Appendix D). The generic evolution of this spectrum under force-free conditions is complex. 
However, it appears to be generic that the evolution of the electromagnetic field will lead to the formation 
of expanding surfaces on which E — > B (cf. Heyl & Hernquist 1998, also Appendix F). (Recall that there 
is no mathematical limitation on B 2 — E 2 changing sign under strict force-free conditions.) In practice, the 
particles that are present are subject to rapid acceleration through E x B drift and this is followed by 7-ray 
emission through synchrotron radiation and inverse Compton scattering, followed by pair production. This 
process will continue until there is enough inertia and or pressure starts to contribute to the equation of 
motion so that relativistic magnetohydrodynamics becomes necessary to describe the behavior of the plasma. 

In the case of polar emission (see below) a possible mechanism of particle acceleration is by inductive 
electric fields developing during electromagnetic collapse near the axis of the flow (Trubnikov 2002). This 
type of acceleration resembles laboratory Z-pinches (e.g. Rudakov & Sudan 1997). This mechanism may be 
considered as reconnection at an O-point. 

In addition to the acceleration mechanisms which are based on known non-relativistic schemes, it is feasi- 
ble that acceleration in relativistic, strongly magnetized plasma may proceed through mechanisms that do not 
have non-relativistic or fluid analogues. Examples of this type of acceleration include kinetic electromagnetic- 
type instabilities of the shell surface currents, as proposed by Smolsky & Usov (1996) and in somewhat 
different form by Liang et al. (2003). Since kinetic properties of strongly magnetized relativistic plasmas 
are only beginning, it is hard to predict acceleration efficiency and particle spectra. Numerical studies in 
the coming years will be most important here. 

In these case dissipation regions will constitute expanding volumes within which there will be a local 
equipartition. These volumes will move with relativistic speeds in the average outflow expansion frame 
defined by the average contact discontinuity and this will influence the observed variability properties. Under 
a broad range of initial conditions, we expect that the density of pairs will increase and the temperature 
will fall until pair production balances annihilation. Given the tenuous, optically thin environment that is 
envisaged, the equilibrium "temperature" will be roughly 10-20keV although a pronounced nonthermal tail 
should also be present in the pair distribution function. The gamma ray spectrum may exhibit a feature 
at this energy which will boosted to roughly ~ 1 MeV. A more careful calculation is needed to see if the 
spectral break sometimes seen at around ~ 200 keV can be reproduced (see also Pe'er & Waxman 2003). 

In this paper we do not commit to a particular acceleration mechanism, but give a necessary qualitative 
description of the possible location of acceleration regions in the frame work of our model and outline possible 
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ways in which magnetic dissipation and particle acceleration may proceed. A more detailed discussion of 
the acceleration and emission physics is deferred to a later paper. We envision two possible locations of the 
emission: (i) along the poles, 9 < 1/r, which we associate with short bursts and (ii) in the of the magnetic 
shell, 9 > 1/r, which we associate with long GRBs. We consider these two possibilities in turn. 

9.2. Polar Emission 

The dynamics of the flow near the polar regions is likely to be unstable to magnetic pinching, giving 
dynamics similar to imploding Z-pinches (e.g. Rudakov & Sudan 1997) where explosive X-ray bursts as well 
as electron and ion acceleration are often observed. 14 In the proposed model the outgoing wave consists 
of currents strongly concentrated near the axis. Such concentration of poloidal currents is likely to become 
unstable due to the development of dynamic and resistive instabilities in the strong current region. The 
fastest growing modes are likely to be pinch, kink and helical modes (see Fig. 7). 

At large currents the plasmas in Z-pinches and TOKAMAKs are known to form thin filaments stretching 
both along and across the direction of the current (Trubnikov 2002). The typical growth rate of filamentation 
instabilities is ~ uj p Vd/c where ui p is the plasma frequency and Vd is the drift velocity of the current particles 
Molvig (1975). Thus, as the current is compressed, at some stage it will break into filaments. Secondary 
filaments will also be compressed and will undergo further filamentation. The development of instabilities 
will result in disruption of the current, during which large DC electric fields are created near the O-type 
point. (The dynamics of current disruption is considered in more detail in Appendix E). This electric field 
will accelerate run away particles (both electrons, positrons and ions). Thus, the magnetic energy will finally 
be radiated as electron cyclotron emission and will also be transported away by cosmic rays. The current 
disruption at each filament will produce a burst of emission, giving a complicated pulse profile for prompt 
GRB emission. 

Initially, the size of the unstable region is of the order of the core radius ro (see Eqns. (30-31)). As 
the instabilities develop the magnetic energy is dissipated, which leads to a loss of magnetic pressure, so 
that progressively larger scale modes become unstable (this is similar to an outgoing rarefaction wave). A 
typical angular size for an unstable current will be of the order of 1 /T since in relativistic flows lateral (in 9 
direction) forces (e.g., magnetic hoop stresses) are strongly suppressed by a factor 1/r 2 (e.g. Bogovalov & 
Tsinganos 1999), so that globally (on scales 9 > 1/r) the flow cannot adjust to the change of equilibrium on 
a dynamical time scale - the time scale for establishing an equilibrium in the lateral direction is much larger 
than the dynamical times scale, by a factor ~ T. The fastest evolution of the flow will take place near the 
axis, within 9 ~ 1/r, where geometric effects may "beat" the relativistic suppression and produce changes 
in the field configuration on shorter time scales. 

9.3. Shell Emission 

The expanding shell will be preceded by a surface (Chapman-Ferraro) current emanating from the poles 
and flowing to the equator (or vice versa). It will be followed, at a distance of the order of the shell thickness 
by a reverse current. There are several possible instabilities associated with this configuration. (Since the 



14 The idea of Z-pinch collapse has already been suggested in application to GRBs (Trubnikov et al. 1995), but the model 
lacked important astrophysical ingredients. 
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shell decelerates on average, it is stable to the Kruskal-Schwarzschild instability). 

Qualitatively, the interface between two magnetized media is known to be unstable in the case of the 
Earth's magnctosphere. Approximately 10% of the incoming energy flux is dissipated at the day side of 
the magnctosphere (e.g. Cowley 1982). In the case of relativistic flows, it is feasible that a similar or 
even larger fraction of the incoming energy flux may be dissipated. Dissipation may be initiated by several 
types of instabilities. First, the surface Chapman- Ferraro current becomes unstable due to development 
of microscopic instabilities on the scales of tens to hundreds of ion Larmor radii (Smolsky & Usov 1996; 
Liang et al. 2003). As a result, the surface current breaks into strongly non-linear fluctuating current and 
charge layers. These current substructures become strongly dissipative and accelerate particles. Secondly, if 
the source varies considerably during its activity the surface current may become unstable due to impulsive 
Kruskal-Schwarzschild instability (see appendix D). The impulsive Kruskal-Schwarzschild instability (IKS) 
develops if the surface separating magnetic field and matter is accelerated by an electromagnetic pulse. It 
is similar to Richtmyer-Meshkov instability in fluid dynamics. At the nonlinear stage, IKS instability will 
lead to formation of small scale structures that will become resistive. During the prompt emission phase, 
the presence of resistive surface currents may lead to an absorption of a large fraction of the incoming fast 
mode energy flux. This will result in dissipation of magnetic energy and particle acceleration. 

An attractive feature of the shell model is that it predicts 7-ray emission over large solid angles and few 
orphan afterglows. The strength of the burst may also depend on the efficiency of radiation generation, which 
may be a strong function of polar angle, but the kinetic total energy (inferred, e.g. , from early afterglows) 
should remain approximately constant. Burst viewed from a small angle (as may be inferred from achromatic 
breaks in the afterglow emission) should be seen to larger distances although a large burst to burst dispersion 
should be expected. In addition, if emission is confined to a narrow region near the surface of the shell, it 
will be more variable than in the case of a filled shell, since there is no "radial averaging" over emitters. 
Simulations of shell emission are also underway which will clarify some of these points. 

9.4. Variability 

One of the key observational properties that every model of GRBs should address is the short times 
scales of variability, as short as ~ 10 ms. One merit of the electromagnetic model is that it allows the 
GRB to originate at a much larger radius, up to <~ 10 16 cm than in the standard, baryonic shock model 
<~ 10 13 cm. This stems from the fast that, in an electromagnetically dominated medium, the emission regions 
(e.g. , large amplitude waves likely to occur in an intermittent turbulence spectrum) may be moving with 
relativistic velocities in the bulk frame. In this section we consider the statistical properties of radiation 
emitted by an ensemble of the emitters confined to a narrow spherical shell and subject to relativistic bulk 
and random motion. We show that relativistic internal motion of " fundamental emitters" can account for 
highly intermittent GRB profile with smaller bulk Lorentz factors or larger radii of emission than is usually 
inferred. 

Consider a shell of thickness AR <C R (in the lab frame) moving relativistically with a bulk Lorentz 
factor r (see Fig. 10). Assume that the shell consists of randomly distributed emitters which move randomly 
with respect to the shell rest frame with a typical Lorentz factor jt- If K is the angle (measured in the rest 
frame of the shell) between the radial direction and emitter's velocity and 1 < 7' < -fx is the Lorentz factor 
of the particle then in the shell frame, emission is beamed into a cone with opening angle A8' r ~ 1/7'. In 
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the observer's frame the emitter's Lorentz factor is 

7 = 7'r(l + ?/Vcos6>;) ~7T(l + cos6»;) (112) 



where v' — y/i — I/7 2 and V = yj\ — 1/r 2 , and angle between the radial direction and the particle's velocity 
in the lab frame is 

tanfl r = - 

r(v + i;cos6';) r v ' 

where the approximation assumes V <~ u ~ 1. From Eq. (113) it follows that as long as u < V, the maximum 
angle with respect to the radial direction that an emitter may have in order to be seen by an observer is 
~ 1/r, regardless of the velocity of the emitter. At the same time, the emission cone in the lab frame is 
A9 r ~ I/7 = l/7rr(l + cos 0' r ). The addition of internal relativistic motion makes the emission beams 
much narrower but the visible emitting volume increases only by a factor of 4, when compared with cold 
emitters. For relativistic internal motion, the number of emitters seen at a given time will decrease as I/77- 
At the same time the flux from each emitter will increase as 7^ (Two factors of 7t coming from the fact that 
all the emission is confined within a narrow beam 1/7^ and two factors of 7^ coming from the relativistic 
contraction of a pulse). Thus, the "thermal"' spread in the motion of emitters can drastically change the 
estimates of the Lorentz factors. The Lorentz factor inferred from variability is, in fact, a product of the bulk 
and "thermal" Lorentz factors. Since for a given T the number of emitters seen at a given time decreases 
oc 1/7^, modest values 7t < T suffice to produce large variations. 



9.5. Causal structure of electromagnetic outflows 

Electromagnetic outflows have a different casual structure from hydrodynamic flows. Initial evolution 
of both types of flow is similar: close to the central source both electromagnetic and hydrodynamic flows 
are subsonic and fully causally connected. Later, after breakout, both types of flows are linearly accelerated 
by magnetic and/or pressure forces r <~ r. Hydrodynamic flows become causally disconnected after passing 
through a sonic point where T = y/3/2. An electromagnetic flow can pass through the fast magnetosonic 
point provided that initially <r > T 2 , and will become causally disconnected over small polar angles AO ~ 

i/r. 

After the photosphere, when magnetization parameter a increases by many orders of magnitude, causal 
behavior of the types of flow is drastically different. Hydrodynamic flows virtually do not establish a causal 
contact over angles large than 1/r. Contrary to that, electromagnetically-dominated flows quickly re- 
establish causal contact. Lyutikov et al. (2003) showed that sub-Alfvenic ejecta re-establishes a causal 
contact over the visible patch of 1/r in just one dynamical time scale (after doubling in radius). In fact, 
causal contact may be established over whole expanding shell after a time t c <~ t s T 2 . At this point the shell 
is fully causally-connected . 

Causal contact over large angles of the shell is established by waves propagating almost backward in 
the shell frame. Here lies a principal difference between strongly magnetized flows on the one hand and 
weakly and non-magnetized flows on the other. Strongly magnetized flows can be subsonic with respect to 
the velocity of the signals in the flow, while weakly magnetized flows are always supersonic. Thus, in case 
of weakly magnetized flows waves emitted backward in the shell frame are "advected" with the flow and 
cannot reach to large angles. Strongly magnetized flows are subsonic, so that waves can outrun the flow and 
reach to large angles. In doing so they leave the shell (which has a thickness either ct s or R/T 2 ), but this 
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does not present a problem since the fast magnetosonic waves become light waves if they propagate into a 
low density medium. 

The causal structure of electromagnetic shells plays an important role in maintaining the coherence of 
large scale magnetic fields. If a surface of relativistically expanding magnetized shell is perturbed at some 
radius, it can quickly propagate information (e.g. magnetic pressure) over the visible angle 1/r, so that the 
shell can have quasi-homogeneous properties despite possible inhomogeneities in the circumstcllar medium 
and in ejecta. In hydrodynamic, a <C 1, or hydromagnetic, a <~ 1 (e.g. Spruit et al. 2000), models only 
under strict homogeneity of the surrounding medium and of the ejecta will the two causally disconnected 
parts of the flow have similar properties. 

Thus, the electromagnetic model provides a solution to the puzzle of how to launch a blast wave that 
extends over an angular scale >> and where the individual parts are out of causal contact. In the 
electromagnetic model, the energy is transferred to the blast wave by a magnetic shell which is causally 
connected at the end of the coasting phase. 



As the magnetic shell expands, its energy is gradually transfered to the preceding forward shock wave. 
In a constant density medium this transfer occurs at the end of the coasting phase. At later times, most of 
the energy released by the central, spinning, magnetic rotator is carried by the shocked interstellar medium. 
The structure of the blast wave is well approximated by the self-similar adiabatic solution (Blandford & 
McKee 1976, 1977), but with angle-dependent expansion. If the total energy released per unit solid angle is 
En ~ L n t s , then (Blandford & McKee 1976) 



As long as the blast wave remains relativistic the energy in the magnetic shell decreases slowly (loga- 
rithmically), so that the relativistic blast wave stage is coexistent with the self-similar stage of the magnetic 
shell as long as the shock remains relativistic, r < tnr = (-Ec/pc 2 ) 1 / 3 <~ 10 18 cm. The afterglow phase 
usually becomes unobservably faint after the expansion speed becomes mildly relativistic. At that point the 
forward shock finally detaches from the magnetic shell and will expand non-relativistically and become more 
spherical with time, resembling a normal supernova remnant. 

Relativistic particles are accelerated in the blast wave producing the observed afterglow in a manner 
which is essentially similar to what is proposed for fluid models. There are, however, some differences. 
First, the contact discontinuity itself may be an important source of magnetic flux through IKS instability 
(appendix D). Then the afterglow may result from a mixture of relativistic particles, derived from the 
shock front with magnetic field derived from the shell, much like what seems to happen in regular supernova 
remnants. 



10. Afterglows 




(114) 



If, in addition, the energy release is as argued above (see Eq. (110)), then 




(115) 



The development of the IKS instability will lead to distortion of the CD. The perturbations propagate 
mostly orthogonally to the magnetic field. In the linear regime, this will lead to "braiding" of the field lines, 



-44- 



so that large scale ordered field remains mostly unaffected. In the nonlinear regime a degree of ordering of 
the magnetic held may be preserved as well. Later, at the afterglow phase, the IKS instability will lead to 
mixing of the ejecta with circumstellar material. This requires that the source remains active for much longer 
than 100 sec. There are indications, e.g. in case of GRB030329, that this indeed may be the case. Since 
during the development of the IKS instability held lines will mostly braid keeping large scale field ordered, 
an ordered component of the field will be introduced to the shocked circumstellar material and may explain 
the polarization of afterglows. 

In addition to dynamic instabilities, the CD may also be unstable to resistive instabilities like the tearing 
mode. The magnetic field from the shocked circumstellar material piles up on the CD, so that however small 
the field is in the the bulk of the flow, it plays an important part in a thin boundary layer near the CD 
(Lyutikov 2002a). Thus the CD becomes a rotational discontinuity and should be susceptible to resistive 
instabilities, similar to the case of Earth magnetopause (Galeev et al. 1986). The growth of tearing mode 
occurs on wave modes propagating orthogonally to the magnetic field. As a result "percolated" magnetic 
filaments form that connect outside and inside plasma. 

The second difference in the afterglow dynamics between the electromagnetic 1 and fireball models is 
that the form of the forward shock is a definite function, determined by the current distribution in the shell. 
This is also the principal difference between relativistic and non-relativistic (Sedov) blast waves. In the non- 
relativistic case explosion quickly "forgets" the details of the central source and dynamics is determined only 
by the total energy release, while in the strongly relativistic regime the forward shock also carries information 
about the angular energy distribution of the source. The form of the blast wave would reflect both the form of 
the driver (magnetic shell) and subsequent evolution of the shock (Kompancets 1960; Shapiro 1979; Granot 
et al. 2003). 15 In the strongly relativistic regime, the form of the forward shock is mostly determined by 
the form of the electromagnetic driver, so that the motion of the shock is mostly ballistic (Shapiro 1979), 
with little sideways expansion (so the in the rest frame v' <C c, not <~ c/y/3 as is commonly assumed); this 
is confirmed by recent semi-analytical (Granot et al. 2003) and fully relativistic hydrodynamic calculations 
Cannizzo et al. (2003). 

Thus, there is little lateral evolution between the forward shock and the contact discontinuity. If we 
continue to use our simple model, we find that the afterglow expansion varies most rapidly and remains 
relativistic for longer times closer to the symmetry axis. In particular, if the outgoing waves correspond 
to the outflow carrying axial current (section 8) the energy carried by the forward shock wave will scale 
as Eq ~ 1/ sin 2 9 and the Lorentz factor of the forward shock T s ~ To/ sin 9. (Note that at the afterglow 
stage the Lorentz factor of the forward shock depends only on the energy that the source released per 
unit solid angle Eq and different from the Lorentz factor of the CD Tcd)- This type of shock has been 
named "structured jet" (or universal jet) (Rossi et al. 2002a), though in our model there is no proper "jet", 
but simply a non-spherical outflow. As Eq, oc sin -2 9 the energy contained in each octave of 9 is roughly 
constant, so that the inferred explosion energy with our simple model will be roughly independent of 9 and 
characteristic of the total energy. This, combined with the assumption that the total energy of GRBs is 
related to kinetic energy of a critically rotating relativistic stellar object, explains a narrow range of inferred 
GRB energies (Frail et al. 2001). Most importantly, the resulting non-spherical blast wave emits in all 



15 Non-spherical shocks and outflows have been extensively studied in the non-relativistic hydrodynamics (Kompancets 1960; 
Laumbach & Probsteinei al. 1969) and in application to the propagation of AGN jets {e.g. Wiita 1978; Norman et al. 1982). 
Shapiro (1979) generalized the work of Kompancets (1960); Laumbach & Probsteinei al. (1969) to non-spherical relativistic 
shock waves. These works are virtually neglected in GRB research. As a result incorrect assumptions {e.g. , "gramophone-type" 
profiles) about the dynamics of the non-spherical shocks were made. 
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directions. However, the intensity of this emission is strongest along the poles. This means that the most 
intense bursts and afterglows in a flux-limited sample will be seen pole-on and can exhibit achromatic breaks 
when T ~ which might be mistaken for jets. The observational appearance of such "structured jets" 
have been investigates by several collaborations (Rossi et al. 2002a; Granot et al. 2003; Perna et al. 2003). 
The overall conclusion is that "structured jets" provide the best fit to afterglow emission properties (Granot 
et al. 2003) and reproduce well luminosity function of GRBs (Perna et al. 2003) including a possibility of 
unified description of GRBs and X-ray flashes (the only X-ray flashes with measured redshift had a total 
energy comparable to classical GRBs (Soderberg et al. 2003)). 

In conclusion, we argue that observational appearance of GRB afterglows depends mostly on two pa- 
rameters: (i) explosion energy (more precisely, on the ratio on the explosion energy to circumstellar density) 
and, most importantly, (ii) the viewing angle that the progenitor's axis is making with the line of sight. This 
possibility, that all GRBs are virtually the same but viewed at different angles resembles unification scheme 
of AGNs. 



11. Observational implications 

Prompt and afterglow polarization. Recent observations by the RHESSI satellite have been in- 
terpreted as evidence for large polarization of the prompt 7-ray emission of GRBs (Coburn & Boggs 2003). 
This has bee contested by (Rutledge & Fox 2003). If the large polarization is confirmed, it could be most 
naturally explained as a synchrotron radiation from a large scale magnetic field. In order to produce very 
high polarization the coherence scale of the field should be larger than the size of the visible emitting region. 
This can be achieved naturally achieved if the flow is magnetically-dominated (Lyutikov et al. 2003) . 

Several natural correlations between the prompt GRB polarization and other parameters follow from 
the model and can be tested with future observations. The maximum amount of polarization in our model 
is related to the spectrum of emitting particles, being higher for softer spectra (Lyutikov et al. 2003). This 
points to a possible correlation between the amount of polarization and hardness of the spectrum. 

Large scale field structure in the ejccta emission may also be related to polarization of afterglows if fields 
from the magnetic shell are mixed in with the shocked circumstellar material. The fact that the magnetized 
boundary becomes "leaky" and both plasma and magnetic field are transported across it has been amply 
demonstrated by decades of space experiments (e.g. Cowley 1982). The transport occurs either due to 
microscopic resistive instabilities of the surface current (similar to the so called flux transfer events at the 
day side of Earth magnetosphere and instabilities at the heliopause (Fahr et al. 1986)) or due to dynamic 
(e.g. RT) instabilities. In the case of Earth and Solar heliopause approximately 10% of the incoming plasma is 
transfered through the CD. In case of GRB the CD may also be unstable due to dynamical impulsive Kruskal- 
Schwarzschild instability (appendix D). Since both resistive (e.g. tearing mode) and dynamics instabilities 
are expected to develop mostly on the waves propagating orthogonally to the large scale magnetic field, a 
degree of field ordering will be preserved. As a result, a comparably large fractional polarization may be 
observed in afterglows as well. In this case, since the preferred direction of polarization is always aligned 
with the flow axis, the position angle should not change through the afterglow (if polarization is observed both 
in prompt and afterglow emission the position angle should be the same). Also, polarization should be most 
independent of the "jet break" moment. This is in a stark contrast with the jet model, in which polarization 
is seen only near the "jet break" times and the position angle is predicted to experience a flip during the 
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"jet break" (Ghisellini & Lazzati 1999; Sari 1999) 16 . Current polarization data are not completely decisive. 
In most cases the position angle remains constant (Gorosabel et al. 2003; Covino et al. 2003a, b; Barth et 
al. 2003; Bersier et al. 2003), while the amount of polarization does not show any correlation with the 
"jet break" . This is consistent with the presence of large scale ordered magnetic fields in the afterglows. 
But there are also exceptions Rol et al. 2003; Grcincr et al. 2003, when the position angle does show some 
variations. (A model of Rossi et al. (2002a) of structured jets also predicts constant position angle, but since 
no large scale magnetic field is assumed the polarization features are still related to the jet break times). 

Structured jet. As discussed in Section 10 our model gives a theoretical foundation for the "structured 
jet" profile of the external shock. We would like to stress once again that there is no proper jet in our model, 
but a non-spherical outflow, so that a break in the light curve is seen when T6 b ~ 1- 

XRF flashes. Another testable prediction of the model is that we should observe much more numerous 
X-ray flashes (XRFs), which may be coming "from the sides" of the expanding shell where the flow is less 
energetic and the Lorentz boosting is weaker. In addition, the total bolometric energy inferred for XRFs 
(presumably from observations of early afterglows before radiative losses become important) should be 
comparable to the total bolometric energy of 7-ray bursts. (The total X-ray fluences may be quite different 
since this would include a correction for the unknown radiative efficiency as a function of the polar angle). 
At present, a single X-ray flash with measured redshift suggests had an energy similar to those associated 
with GRBs (Soderberg et al. 2003). Generally, the distributions of parameters of XRFs should continuously 
match those of GRBs. 

Short-long dichotomy. Our model can also explain a short versus long dichotomy in GRBs. We 
associate short bursts with the instabilities of polar currents while long with instabilities of shell currents. 
Our analytical solutions for electromagnetic fields diverge on the axis; in reality the axial current density 
will have finite angular width. If an observer is located within this small angle he would see a burst that 
would be shorter since the Lorentz factors of the outflow arc higher closer to the axis. 

Hard-soft evolution and E pea k — \[L correlation. Magnetically-dominated flows may also explain 
the observed correlations between various GRB parameters, though the lack of a testable particle acceleration 
model makes the arguments below only suggestive. For example the trend of GRB spectra to evolve from 
hard to soft during a pulse is explained as a synchrotron radiation in an expanding flow with magnetic field 
decreasing with radius B oc \/~L/r (later in a pulse emission is produced further out where magnetic field 
is weaker, so that the peak energy will be lower; this is similar to " radius-to- frequency mapping" in radio 
pulsars and AGNs). A correlation between peak energy and total luminosity, E pea k ~ \fL (Lloyd- Ronning 
& Ramirez-Ruiz 2002) follows from the assumption of a fixed typical emission radii and fixed minimum 
particle energy. Both of these correlations are independent of the bulk Lorentz factor, but depend on the 
lower energy cut-off in the spectrum of accelerated particles. 

No reverse shock. Another consequence of the model is that no emission from the reverse shocks is 
expected in electromagnetic models. When plasma magnetization becomes of the order of unity relativistic 
MHD shock conditions are modified: shocks becomes less dissipative; in the force-free limit shocks cannot 
exist at all. Early optical flashes, which are conventionally associated with reverse shock {e.g. Mcszaros & 
Rees 1997; Kobayashi & Zhang 2003; Kumar & Panaitescu 2003), should have a different origin. A possible 
explanation is the formation of a radiative precursor that pre-accelerates the medium, loads it with e+- pairs 



16 Note that these works did not take into account kinematic depolarization (e.g. Lyutikov 2003a) which may give an error 
of the order of unity. 
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and produces soft emission (Beloborodov 2002; Thompson & Madau 2000) 

12. Discussion 

In this paper we have explored the "electromagnetic hypothesis" for ultra-relativistic outflows, namely 
that they are essentially electromagnetic phenomena which are driven by the energy released by spinning 
black holes or neutron stars and that this electromagnetic behavior continues into the source region even 
when the flows become non-relativistic. In this model the energy to power the GRBs comes eventually from 
the rotational energy of the progenitor. It is first converted into magnetic energy by the dynamo action of 
the unipolar inductor, propagated in the form of Poynting flux-dominated flow and then dissipated at large 
distances from the sources. We have taken an extreme view that the flows become essentially force-free and 
explored the consequences of this assumption. Real GRB outflows must contain both baryonic matter and 
approach near-equipartition in the emission regions. 

This model envisions GRB outflows as an electric circuit in which the central source acts as a power- 
supply generating a very simple current flow - along the axis, the surface of the shell and the equator. 
In practice, it is the detailed electrodynamics in the vicinity of the outgoing light surface that fixes the 
poloidal magnetic field and electrical current distributions. We can therefore change these (still, of course, 
maintaining the force-free condition) and solve for a new evolution of the magnetic shell and blast wave. 
A broader distribution of currents will generally produce a less pronounced expansion along the axis and 
change somewhat the statistics of the observed afterglows. Alternatively, the central source can collimatc the 
flows to even more narrow (cylindrical) expansion. In the solution above, we ignored the poloidal magnetic 
field and, consequently, the angular momentum. These can be reinstated perturbatively into the solution. 
Their influence wanes with increasing radius. Other ways to obtain different solutions include changing the 
temporal variation of the source from the one satisfying self-similar expansion to a more general variation. 
As the individual parts of the blast wave expand essentially independently, when ultra-relativistic, there are 
no new issues of principle to address in solving these problems. When the blast wave becomes non-relativistic 
the interior gas will be roughly isobaric and the shell will become more spherical with time. 

The most striking implications of the electromagnetic hypothesis are that particle acceleration in the 
sources is due to electromagnetic turbulence rather than shocks and that the outflows are cold, electromag- 
netically dominated flows, with very few baryons at least until they become strongly dissipative. The major 
drawback of the current model is the lack of the detail model of energy dissipation and particle acceleration in 
relativistic magnetically-dominated medium. This is a difficult problem. The problem with magnetic dissi- 
pation may be exemplified by the better studied way of dissipation of magnetic energy, through reconnection. 
Reconnection physics is highly uncertain: it depends crucially on the kinetic and geometric properties of 
the plasma, which is very hard to test observationally. Often, various instabilities (based on inertial or cy- 
clotron effects, on ions or electrons) often seem to account equally well (with astrophysical accuracy) for the 
observed phenomena. This is an important uncertainty since the principle scale of the reconnection region 
(related, for example, to electron or ion skin depth, Larmor radius or magnetic Debye radius) is crucial in 
determining the rate of reconnection (e.g. Birk et al. 2000). This situation may be contrasted with the 
shock acceleration schemes, where a qualitatively correct result for the spectrum of accelerated particles can 
be obtained from simple macroscopic considerations. 

An important implication of the electromagnetic model is that supernova explosions may be magnet- 
ically driven as well (Leblanc & Wilson 1971; Bisnovatyi-Kogan 1971). A number of research groups are 
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investigating this possibility, mostly through numerical modelling (Mizuno et al. 2003; Proga et al. 2003). 

Another uncertainty in the model is how baryon-free the flow can become. At early stages large fluxes 
of low energy neutrinos are expected. They will drive the material up the field lines into the magnetosphere, 
while gravity will try to pull the matter back into the disk. We require that the resulting outflow be extremely 
clear of baryons (~ 10~ 9 in energy flux); numerical simulations are required to see if this is possible {e.g. 
Thompson et al. 2001). 

If GRBs are electromagnetically-driven they should not be associated with strong, high energy neutrino 
sources. A gravitational wave signal is expected in some, though not all, source models and would be strongly 
diagnostic if ever detected. 

Electromagnetic models are also consistent with being possible sites of UHECR acceleration (Vietri 1995; 
Waxman 1995; Milgrom & Usov 1995; Wick et al. 2003) if the corresponding energy balance of UHECR and 
GRBs agree (Waxman 2002; Vietri et al. 2003) and if assuming that the GZK cut-off is indeed observed in the 
spectrum of UHECRs (Bahcall & Waxman 2002). Primarily, the total available electric potential in GRBs, 
<~ 10 22 V, is sufficient for acceleration of UHECRs. UHECR energies are mostly likely radiation-limited and 
not acceleration-limited, but, at late stages of GRB expansion, synchrotron losses become sufficiently small. 

In addition to Fermi acceleration at the forward shock, electromagnetically-dominated flows provide 
several possible acceleration schemes in their own right. (Note that, in order to explain UHECRs, shock 
acceleration models have to assume that turbulent EMF reaches approximate cquipartition behind relativistic 
shocks (e.g. Dcrmer & Humi 2001); this requirement makes these models akin to electromagnetic models.) 
There is a large inductive electric field Eg which creates a large potential difference (4) between the polar 
and equatorial flux surfaces. Thus, in order for a particle to gain the total energy it should traverse from 
pole to equator (or a considerable part of this path). This should be done on expansion time scale in order 
to avoid adiabatic loss; this is possible at the end of the relativistic stage. (Note that, at late stages, most 
of the energy is in the forward shock the total pole-to-equator electric potential drop remains approximately 
constant as long as the expansion is relativistic and starts to decrease only at the non-relativistic stage.) 
This crossing of the potential surfaces may be done either kinetically, e.g. due to a drift, or resistively, e.g. 
at O-type point (e.g. Trubnikov 2002; Vasyliunas 1980). 

In this paper we have discussed the basic principles that may be important in electromagnetically- 
dominated outflows. At the present stage the model is definitely less developed than the hydrodynamics 
fireball model and a lot needs to be done to vindicate it, but the possible fundamental advantages of electro- 
magnetic models outlined here make it, in our opinion, a very promising approach. From a more theoretical 
perspective, there is much to be learned about the properties of force-free electromagnetic fields, especially 
their dissipation. The possible relationship of the GRB fluctuation power spectrum to an underlying turbu- 
lence spectrum is especially tantalizing. Undoubtedly, numerical simulations will be crucial as the problem 
is essentially three-dimensional. Force-free electromagnctism is easier to study than relativistic MHD and 
may well be a very good approximation in many of these sources. 

If, as we have argued, magnetic fields do play an important dynamical role in GRBs, there is the between 
supersonic MHD outflows a < T 2 (e.g. Vlahakis & Konigl 2003a; Fendt & Ouyed 2003) and the subsonic 
force-free model, a > T 2 , considered in this paper. As we have discussed these two alternatives have very 
different properties coming from different causal structures of the flow. Obscrvationally, it will be hard to 
distinguish between the two since even in the strongly magnetized case we expect that in the emission region 
the fields approach cquipartition (since then plasma emits most efficiently). 
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Finally, we give a concise list of the observational properties of the GRBs which are explained or foreseen 
in the electromagnetic model. GRBs properties should most depend on the observer angle with respect to 
the jet axis. Electromagnetic model produces "structured jet" with energy Eq oc sin~ 2 9; there is no problem 
with "orphan afterglow" since GRBs are produced over large solid angle; X-ray flashes are interpreted as 
GRBs seen "from the side" , but their total energetics should be comparable to proper GRBs; overall energy 
of GRBs is expected to have a small scatter if it is related to a critically rotating relativistic object; the 
model can qualitatively reproduce hard-to-soft spectral evolution as a synchrotron emission in ever decreasing 
magnetic field B oc VL/r (L is luminosity, r is emission radius), akin to "radius-to-frequency mapping" in 
radio pulsars and AGNs; similarly, the correlation E pea k ~ \[L is also a natural consequence. Finally, high 
polarization of prompt emission may also be produced: it should correlate with the spectral index; mixing 
between circumstellar material and ejecta will result in the same position angles of the prompt emission and 
afterglow which should be constant over time; fractional polarization should be mostly independent of the 
"jet break" time, but may show variations due to turbulent mixing. 

13. Application to other sources 

Many of the principles described in this work may be applied to other astrophysical sources like pulsars, 
(micro)quasars and AGNs (Blandford 2002). Despite the apparent differences (inner boundary condition for 
neutron stars and black holes, time (non-)stationarity, optical thickness) it is plausible that all these sources 
produce ultra-relativistic magnetically-dominated outflows with low baryon density. Energy, transported 
primarily by magnetic fields, is dissipated far away from the source due to the development of MHD current 
instabilities. Particles are accelerated in localized current sheets by DC electric fields and/or electromagnetic 
turbulence producing bright knots (in AGNs) and a variety of bright spots in pulsar jet, best observed in 
the Crab. In situ acceleration is required in both AGNs and pulsars since the life-time of emitting particles 
is shorter than the dynamical time. Qualitatively, magnetically-dominated jets are expected to be more 
stable than their fluid counterparts (Bcskin & Pariev 1993). In the case of AGN, recent observation of TeV 
emission from blazers with a very short time scale variability (Krennrich et al. 2002) favors magnetically- 
dominated jets {e.g. Lyutikov 2002b) since the two leading models of jet composition - leptonic and ionic - 
both have difficulties in explaining the fast variability (problem for ion-dominated jets, since cyclotron times 
are very long for ions) and total energy content (since leptons suffer strong radiative losses at the source 
due to radiation drag). The flow evolution in all these systems may proceed in a similar wave. Jets are 
launched in the vicinity of a central object (a BH-disk system or a neutron star) where the flow is subsonic 
and in a dynamical balance. As the flow propagates out jets become laterally imbalanced, just as with 
GRBs, the lateral dynamics is suppressed by relativistic kinematics, so that the jet cannot adjust to the 
change of equilibrium on dynamical time scale. After a jet has propagated T 2 dynamical times, the lateral 
dynamics starts to become important. A particularly strong evolution will take place in the core of the 
jet, where magnetic hoop stresses may lead to the collapse of the flow toward the axis. For example, AGN 
jets may originate close to the central black hole with a Schwarzschild radius r s <~ 10 13 cm with a typical 
Lorentz factor T ~ 10 — 30. Then the dynamical evolution and the onset of emission occurs on a scale 
~ r 2 r s ~ 10 15 — 10 16 cm. In pulsars acceleration occurs on a scale of light cylinder, r s ~ 10 10 cm, with a 
Lorentz factor T <~ 100 — 1000, then the dynamical instabilities occur at <~ 10 14 — 10 16 cm. 

Over the past three and a half years that this work has been in progress we benefited from numerous 
discussions with our colleagues. In particular, RB thanks Jerry Ostriker for early discussions of electromag- 
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netic models of gamma ray bursts and ML thanks Chris Thompson for discussion and comments on the 
manuscript. This research has been supported by NSERC grant RGPIN 238487-01 and NASA grant 5-2837 
and 5-12032. 
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In this appendix we consider the far field dynamics of a spherically symmetric magnetized wind. The 



the electro-motive force £ produced by the central source. The central source loses energy in two forms: 
mechanical Lm and electromagnetic Lem luminosities. 



We wish to understand how the parameters of a fully relativistic flow (velocity /?, pressure p, magnetization 
cr) evolve for an arbitrary ratio of both Lem/Lm and p/p (p is the rest- frame mass density). 

The asymptotic evolution of the flow is determined by conserved quantities which may be chosen as the 
total luminosity L, the mass flux M and the EMF £. Thus, the central source works both as thruster and as 
a dynamo. We follow the flow evolution starting from some small inner radius where we assume that a flow 
with given L, M and £ (in a form of toroidal magnetic field) is generated. Thus, we avoid the important 
question of how the flow is launched. Another important assumption that we make is that the expansion is 
radial. 

The formal treatment of the problem starts with the set of RMHD equations which can be written 
in terms of conservation laws (50) in which we assume that fluid is polytropic with adiabatic index T a : 
w = p + r r ° 1 p. (We assume that the adiabatic index T a is constant for algebraic simplicity.) 

Writing out eqns. (50) in coordinate form and assuming a stationary, radial, spherically symmetric 
outflow with toroidal magnetic field, we find 



A. Dynamics of magnetized pair-loaded flows (minifireballs) 



dynamics of such warm magnetized wind is controlled by three parameters: energy L, mass flux M and 



L — L M + L E m 



(Al) 




(A2) 



la r .[r 2 (( W + 6 2 )/3 2 r 2 + (p + 6 2 /2))] 
-d r [rb[3T] = 



r 




(A3) 



(A4) 



\dr [r 2 p(5T] = 



(A5) 



The above relations can be simplified if by defining 




a 



(A6) 



where £ is the electromotive force. 



This preprint was prepared with the AAS IATjTjX macros v5.0. 
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It is convenient to introduce two other parameters: the magnetization parameter a as the ratio of the 
rest-frame magnetic and particle energy-density and a fast magnetosonic wave phase velocity (3f 

Using the three conserved quantities L, M and £ the evolution equation becomes 
i (r o -l) ((3L-f3TM-£A 

d r r = — i —L r (A8) 



'2 



2/3 ' r r (/3L(/? 2 + 1 - r o ) + (T„ - 1)0TM - (2 - T a )£ 

Eliminating £ in favor of T f we get a particularly transparent form for the evolution of Lorentz factor 

( r2 - r / 2 ) g r - 2pTa (A9) 

Equation (A9) is nozzle- type flow (e.g. (Landau & Lifshits 1975a)); in astrophysical context it is best known 
for Parker's solutions of the solar wind (Parker 1960). The lhs of eq. (A9) contains a familiar critical point 
at the sonic transition T = Tf. The positively defined first term on the rhs describes the evolution of Lorentz 
factors due to pressure effects. In the case of purely radial expansion the magnetic gradient forces are exactly 
balanced by the hoop stresses, so that magnetic field does not contribute to acceleration. From Eq. (A9) it 
follows that super-fast-magnetosonic flows accelerate while sub-fast-magnetosonic flows decelerate. It also it 
follows that terminal Lorentz factor of the flow is determined by the condition d r T = which implies that 
either p = or = 0. Condition (3^ = can be reached only for subsonic flows with £ = 0; (a steady 
state perfect MHD radial magnetized flow with £ 7^ cannot slowdown to a halt at infinity). Neglecting the 
floo = solution, the terminal velocity of magnetized flow is determined only by the condition p = 0: 

£ 2 

L = r^M + — (A10) 

Poo 

For each set of parameters (L, M and £) there are generally two solutions for the terminal four- 
velocity corresponding to the supersonic and subsonic solutions. Th only exception is the unmagnetized 
branch, £q. In the absence of magnetization the terminal (supersonic) Lorentz factor is uniquely determined 
by Too = L/M. In GRB applications we are mostly interested in the supersonic solutions (the subsonic 
solutions are important for pulsar winds, {e.g. Kennel & Coroniti 1984)). 

For supersonic solutions, the unmagnetized branch with £q is qualitatively different from magnetized 
branch. In the absence of magnetization, the terminal (supersonic) Lorentz factor is determined by Too = 
L/M (in the subsonic regime f3 ~ r~ 2 and p, p = const as with the breeze solution of the solar wind). For 
finite magnetization there are two solution: a larger one corresponding to super-fast-magnetosonic flow and a 
smaller one corresponding to a sub-fast-magnetosonic flow. In particular, for large L/M>1 the supersonic 
terminal proper velocity reaches 

TooPoo ~ (All) 

Thus, for a given ratio L/M the terminal Lorentz factor decreases with increasing £ - magnetic field provides 
an effective inertial loading of the flow. 

For non-zero magnetization the terminal velocity cannot be determined uniquely from a given M and 
£ - this problem is the results of our neglect of details of acceleration (Michel 1969; Goldreich & Julian 
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1970). Yet, there is a trick, initially suggested by (Michel 1969), which allows one to select from a family of 
solutions a "minimum torque" solution (in the terminology of Michel (1969)). A more careful examination 
(Goldreich & Julian 1970; Kennel et al. 1983) proves that indeed, for strongly magnetized flows, when the 
pressure contribution to acceleration is not important, the minimum torque solution is the correct one. 

In case of non-zero magnetization, for a given ratio £ solutions exist only for L/M larger than some criti- 
cs/ 

3/2 



calvalucL/M= (l-(£ 2 /i) 2 / 3 )- 3 / 2 , corresponding to /3oo = (£ 2 /L) 1 / 5 (alternatively, = l/yjl + (M/£ 2 ) 2 / 3 , 

/ . \3/2 

reached at £ 2 1 1 + (M/£ 2 ) 2 ' 3 J ). This is the Michel solution, which we can reformulate now stating that 

for fixed M and £ the minimum energy loss is reached at /3 m in- It is expected that, as the flow become 
stronger magnetized, its terminal evolution asymptotes to the Michel solution. 

The above relations may be expressed in term of terminal magnetization parameter and the terminal 
velocity of the supersonic wind 

L £ 2 

Toe = -. ,, (Toe = — . (A12) 

M (l + a oo )" T^PooMo 

Assumption floo = m i n , then gives 

(l + CSooroo^oo) 2 / 3 )^^^^ (A13) 

Poo 

which in the strongly relativistic limit gives the Michel solution 

Too = (A14) 

We can also relate the terminal magnetization <7oo to the magnetization at the source - more specifically 
to magnetization at the sonic point 07. Using (A12) we find that at any point in the flow the magnetization 
parameter is given by 

u = — (A15) 

/3(l+<7oo)-/3oo<Toc 

To relate the magnetization at infinity to the magnetization at the sonic point we need solve for the velocity 
at the sonic point /3/. Using (A7) and (A15) we find an equation for (3f 

0} = (T. - 1) (l - = r ) + (2 - T a )-^-^ (A16) 

V loo(l + CToc)/ l + cr 



In the absence of mass flux, M = 0, Too = 00, 



f^+yfr „ << ; (Aly) 

Thus, for strongly magnetized flow, (Too ^> 1, T/ = ^/2a oc for T a = 4/3. Using above expressions for the 
sonic velocity we find a relationship between the magnetization at the sonic point and infinity (for supersonic 
flows) : 

if (Too < 1 



Thus, we always have < <Jf, but they remain of the same order of magnitude: the magnetization of the 
flow changes only slightly as the flow propagates away from the launching point to infinity. The reason for 
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constant a in the supersonic regime is that both in the case p > p (linear acceleration stage) and p <C p 
both the plasma and the magnetic field energy densities in the flow change with the same radial dependence 
(~ r~ 4 and <~ r~ 2 correspondingly). 

Generally there are two branches of solutions: supersonic and subsonic. In the limit — > oo (this is 
equivalent to neglecting the mass loss rate of the central source in comparison with the energy loss rate), 
the nozzle equation (A8) simplifies 

(3 lnr " & (/3 2 - r a + 1) - (2 - r a )£*/L (AW) 

which can be integrated to give 

r oc (/3r) (2 - r<,)/(2(r °- 1)) (0 - ^r V(2(IW)) = (/3 _ §j W j- 2 for Y a = 4/3. (A20) 

Thus, the supersonic flow is accelerated by pressure effects as long as p » p, reaching a coasting phase with 
r <~ Too when p < p. while subsonic flows reach a minimum velocity given by 

/?=-^_. (A21) 

For arbitrary flow parameters the evolution equations are integrated numerically (Fig. 4). Given the evolu- 
tion of the flow and the relation for local a (A15) we can find the evolution of the magnetization parameter 
(Fig. 5). 



B. Shocks driven by magnetic explosion 

As the magnetic shell expands into an ambient gas a shock forms ahead of it. The shock would quickly 
go into self-similar regime described by the Blandford-McKee solution (Blandford & McKee (1976)). In this 
appendix we relate the two self-similar solutions - inside the magnetic shell and unmagnetized Blandford- 
McKee solution for the flow between the forward shock and the CD. It is straightforward to generalize these 
results to include external magnetic field using the solutions for the self-similar structure of magnetized blast 
waves derived by Lyutikov (2001). We also neglect the lateral dynamics of the shocked material and assume 
that between the blast wave and the CD (which are separated by ~ -R/r 2 ) the plasma moves radially. This 
is justified as long as the Lorentz factor changes on angular scales larger that 1/r 2 . 

For a self-similar blast wave the contact discontinuity is located at a fixed \ (The self-similar variable 
X in this section is chosen in such a way that \ — 1 011 the- forward shock). Differentiating Eq. (C7) with 
respect to time and using r'(t) = 1 — l/(27 2 ) = 1 — l/(yg) we find 

(l + 2(l + m)y) 2 ^2 

Xcd = t ~ - 

gy (l + 2(l + m) 2 y) 9 

(c.f. Blandford & McKee (1976) Eq. (39)). 

We have to balance the pressures on both sides of the CD. The momentum flux of the shocked material 
in the frame of the CD is 

(4 7 2 - 1) p0 w 4/3T 4 g( X cD)f(xcD)w cxt (B2) 
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where w cxt = n cvit mi is the enthalpy ahead of the shock, and g(xco) and J{xcd) and the Blandford & 
McKee (1976) values on the CD. 

Using Eq. (B2) we find 

g 

W-i = KW cxt , K = -g(xcD)f(XCD) (B3) 

for the enthalpy on the CD. Relation (B3) determines the boundary condition on the CD and, thus, the 
normalization for the solution inside the CD. 

The relations for the g(xco) and J(xcd) are easily found from Blandford & McKee (1976): 

16 + 28 m-3 k m+m 2 2 + n 

/n - 9 h\ 

9(Xcd) 



/a-2b\ " /2(12-3fc-m)Y 

\a + 2b) \ 3(3- k-m) J 

64+3 fc 2 +2 fc (-16 + m)-24 m-m 2 

a + 2b\ 4b ~~ /2(12-3fc-ra) 



f ^={— h ) { 3(3-*-m) ) ™ 

where 

a = 32-9fc-5m, b 2 = 4(l + m)+ (12 ~ 3 4 fc + m) (B5) 

where k is the power law of the density variation p <~ r~ k (see Fig. 8). These relations allow us to sew the 
two solutions - inside the CD and between the CD and forward shock. 



C. Self-similar expansion of magnetic shell: Relativistic MHD approach 

In this appendix we illustrate the similarities of the RMHD and RFF approach re-deriving the self- 
structure of the magnetic shell from RMHD equations. The self-similar approach is useful when the thickness 
of the magnetic shell is much smaller than it radius. For a power-law dependence of Lorcntz factor of the 
shell on time, V cx t~ m / 2 the the fields are concentrated near the CD for 1 < m < 3. In this case we can 
neglect the divergence of the characteristics and apply the hydrodynamic Blandford-McKee approach. 

Writing out Eqns (50) in coordinate form and assuming an azimuthally symmetric outflow with toroidal 
magnetic field, the conservation of energy and momentum, induction equation and mass conservation give 

d t [(w + 6 2 ) 7 2 - (p + & 2 /2)] + \d r [r 2 (w + 6 2 )/3 7 2 ] + — ^—dg [sin0(tu + 6 2 ) 7 2 A] = (CI) 



r 2 L J rsin^ 

d t [(w + b 2 h 2 p] + \d r [r 2 ((w + b 2 )(3W + (p + b 2 /2))] 



l _p \Ji(u.. , ,2^,2 , , 2p+(w + b 2 h 2 X 2 



+ —^d [sin 0(w + 6 2 ) 7 2 A/3] = (C2) 
r sin L J 

d t [(w + b 2 h 2 X] + ±d r [r 2 ( W + fe 2 ) 7 2 A/3] + ^±^!M _ cot B V_111 + 

— [sin 9 ((w + 6 2 ) 7 2 A 2 + (p + 6 2 /2))] = (C3) 
v sin /J 

d t [by] + -d r [rb(3y] + -d [byX] = (C4) 
r r 

d t [pi] + y 2 d r [r 2 p(3y] + -±- e d e [pXy] = (C5) 
where A = vg. 
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Our goal is to find relativistic, self-similar solutions to Eqns (C1-C5) in the limit w, p,p <C b 2 and A <C 1. 
We assume that matter inertia and pressure can be neglect and set w = p = 0. The dynamical equations, 
which are equivalent to force-free equations (40, 43), become 

dt [b 2 j 2 b 2 /2] + \d r [r 2 6 2 /3 7 2 ] + —do [sinflfcVA] = 

dt [bWP] + P (*W + b 2 /2)] + [sin^VA/?] = 

d t [6 2 7 2 A] + \d r [rV 7 2 A/?] + + cot( £ + _J_^ [sintf (& 2 7 2 A 2 + fe 2 /2))] = 

V V ZiT V Sill u 

d t M + -d r [rbpj] + -d e [67 A] = (C6) 
r r 

The strongly relativistic solutions that we are interested in involve expansion of all quantities assuming 
large 7-factors. In addition, assuming that A7 <C 1 and balancing the powers of A in eq. (C3) we conclude 
that in order for the self-similar solutions involving radial and longitudinal expansion to exist, we need 
A oc I/7 2 . 

Following Blandford & McKee (1976) we chose the self-similar variable 

X = 1 + 2(m + 1)£ = [1 + 2(m + l)r 2 ](l - r/t) (C7) 

where £ = (1 — r/i?)F 2 , R = t(l— l/(2(m + l)r 2 )) is the radius of the contact discontinuity and we assumed 
that the Lorentz factor scales with radius as T 2 oc t~ m . We limit ourselves to the strongly relativistic case 
expanding all relations to the first order in 1/r 2 . 

Treating (x,y), where y = T 2 , as new independent variables we find 

d t = -mydy + ((to + l)(2y - x) + l)# x 
d r = -(l + 2(m + l)y)d x 

2yg 

r = t(l- - (C8) 

The boundary conditions on the CD require that Lorentz factor of the magnetic field lines equals the 
Lorentz factor of the CD itself, while pressures on both sides of the CD should be equal. This and the scaling 
A ~ 1/r 2 , allow the following parameterization: 

7 2 = T 2 g( X ) 
b = s/2w[Th(x) 

a = i(x)/r 2 (C9) 

with 5 (1) = h(l) = 1 and T = T(6). 

The equations for the self-similar variables are 17 . 
d In g 1 — m 



d In x 1 + to 



(CIO) 



^General relations for radial motion (I = 0) including effects of nonzero pressure are derived in Lyutikov (2002) 
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d\n X 2(1 + m) 1 ' 

(1 + m)(l - ^ = - sin%1/2 J + (m - 1 - (m + 1) X5 ) 1 (C12) 

Taking into account boundary conditions, the solutions are 

g = X ^, (C13) 
h = x^it^) (C14) 

(compare with 87). 

One particular case is I = - purely radial motion with parameters depending on the polar angle. In 
this case from (CI 2) it follows 

v = r 2 = (Ci5) 

sin 

where To = T 2 (tt/2). This solution corresponds to the fully balanced case (80). 

If the magnetic shell expands into the progenitors wind with density varying as a power law of radius, 
pi oc r~ k the relations (C12) may be easily generalized. For the laterally balanced solution (I = 0) we find 



g dx 4(m+l) 
with solutions 

2(l-m)-fc 
= Y 2(l + m ) 



1 9 In 5 _ 2(1 - to) - k 

(C16) 



g d X ~ 2(m+l) 
1 a In ft 2(1 -to) -A; 



.9 = X 

2(l-m)-fe 

ft = x 4(1+m) (C17) 



It is required for consistency that k < 4 and to > — 1. 



D. Instability of the contact surface: Impulsive Kruskal-Schwarzschild 

In this appendix we consider the dynamical instability of the CD due to a fluctuating luminosity of 
the central source. Conventionally, the Kruskal-Schwarzschild (KS) instability refers to an instability of a 
plasma supported against gravity by magnetic field. If plasma is "below" the magnetic field (so that the 
effective gravity is directed from the magnetic to the plasma phase) then the configuration is stable. This is 
similar to the Ray leigh- Taylor (RT) stability/instability. 

A stable contact discontinuity separating two fluids may become unstable when a shock passes through 
it. In fluid dynamics this is known as Richtmyer-Meshkov (RM) instability. It is independent of the whether 
the contact surface is RT stable or unstable. Physically, small ripples of the CD distort the flow of the 
shocked material and create vorticity which destroys the CD. Mathematically, the RM instability is treated 
as an impulsive acceleration of the CD during the passage of the shock. Similarly, we expect that the KS- 
stable CD will be unstable under influence of an impulsive electromagnetic perturbation propagating in the 
magnetically dominated medium. This instability may be called impulsive Kruskal-Schwarzschild instability, 
(or, equivalently, magnetic Richtmyer-Meshkov instability). 
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The mathematical analysis of the Richtmyer-Meshkov instability is complicated, yet a simple original 
derivation of Richtmyer (1960) gives a good estimate of the growth rates in a wide range of cases (see also 
Inogamov 1999). Below we follow the logic of Richtmyer (1960) in estimating the growth rate of impulsive 
KS instability. A more detailed analysis will be reported elsewhere. 

Consider a pressure balanced CD separating regions of magnetic field and plasma, so that B 2 /8n = p. 
If the perturbation of the CD with a wave number k is propagating orthogonally to the direction for the 
magnetic field the behavior of the displacement £ is determined by 

§-9K = (Dl) 

In the case g > 0, perturbations grow exponentially: this is the KS instability. In the case g < 0, this describes 
oscillatory perturbations propagating orthogonally to magnetic field (akin to surface gravity waves). 

Assume next that the magnetic pressure suddenly increased by SB. This may be due to (fast) magnc- 
tosonic wave pulse propagating towards the boundary. This will launch a shock wave in the plasma, so that 
the CD will acquire some velocity Av. In this case the effective gravity on the CD is g — Av6(t). Integrating 
Eq. (Dl) once we find 

§ = Avk^o (D2) 
at 

where £o is the initial perturbation. Thus perturbations will grow. Growth is initially linear in time (as 
oppose to exponential in case of SM instability). This is the impulsive Kruskal-Schwarzschild (IKS below) 
instability. It's growth rate is 

~ Avk (D3) 

We also note, that IKS, like RM, is not a classical instability: it does not involve "self-amplification", 
so that at the linear stage perturbations grow linearly, not exponentially, with time. The reason is that the 
initial pulse generates velocity fields which evolve dynamically and distorts the CD by inertial motion, so 
that the typical velocity of a perturbation will decrease with time. 

Next we consider the IKS instability in GRBs. The pressure balance in the frame of CD gives (58) 

" =p' = npT 2 (D4) 



where p' is the pressure of the shocked material and k is a constant relating the pressure on the CD to the 
external density. In case of self-similar motion the value of k can be easily found from the Blandford-McKee 
solution (appendix B). 

Magnetic field B(t') at time t' is related to the retarded time luminosity of the source t' ret : 



B{ t ') = -J^lA (D5) 

r V c 

As we have discussed in Section 7, for constant luminosity these equations determine self-similar evolution 
of the magnetic shell. Since Lorentz factor decrease with time the effective gravity in the frame of the CD 
is directed from the magnetic shell ("light fluid") to shocked plasma ("heavy fluid"), so that the system is 
KS stable. 

Next assume that at a retarded time t' ret the luminosity of the source increases by SL so that 

L = L + 5L@{t' ret - t' 0iret ) (D6) 



-65 - 



(O is Hcaviside function). When the fast magnetosonic waves propagating information about the change in 
the source luminosity reach the CD, it will launch a forward shock (or a sound wave) in the already shocked 
circumstellar medium. As a results the velocity of the CD will change. The changes in the velocity may be 
related to the change in luminosity 

SL SB ST r , , , N 

— = 2— = 4— = ASvT 2 /c = A8v'/c (D7) 

Li JD I 

where Sv' is a change in the velocity of the CD in its frame and we assumed that all the changes are small, 
< 1. 

Note, that the electromagnetic pulse associated with the increase of luminosity should not necessarily 
launch a shock. All that is needed is a sharp velocity change of the CD. In the hydrodynamic case such a 
sharp velocity change can be only due to a shock (any subsonic flow would communicate information pre- 
accelerating the CD). In the case of force- free fields, a sharp velocity change may be due to an electromagnetic 
pulse propagating with the speed of light. 

We can now estimate the growth rate of IKS instability in GRBs. From (D3) and (D7) we find the 
growth rate in the frame of the CD: 

JiKS-\ S -^ck (D8) 

Qualitatively, for SL/L ~ 1 

i'iks ~ f (D9) 

Thus, the IKS instability grows on time scales of the order of sound crossing time. 

At the largest scale, r 1 ~ r/T, the instability will be suppressed by the spherical expansion of the flow. 
This will occur when jiks ~ 3r'/r' ~ 3/t' (c.f. the Hubble flow). In the observer's frame the largest 
unstable mode will have a scale somewhat smaller that the "horizon" < r/T. 

In this appendix we briefly discussed the instability of the CD due to source non-stationarity, assuming 
that the source changes its luminosity on the scales much shorter that the dynamic scales of the outflow. 
This is indeed expected since the source may change on a scale of milliseconds. On the other hand, if the 
source luminosity changes on a scale comparable to the outflow time scale, the impulsive KS instability 
becomes a conventional KS (in)stability. The relation between the two is very much similar to the relation 
between RT and RM instabilities and is quantified by the Froude number Fr = Avy^k/g. For RT and KS 
instabilities Fr ~ 1, while for RM and IKS Fr > 1. 



E. Disruption of polar current 

In this Appendix we consider disruption of the axial current due to the development of dynamic insta- 
bilities near the polar axis. As we have argued, the lateral (in the 9 direction) dynamics is suppressed in 
relativistically expanding flows. As a results, at large distances from the source, when the Lorentz factor has 
decreased considerably, the flow will be dynamically unbalanced. In this case the development of dynamic 
instabilities may proceed in an explosive fashion, so that the axial current is disrupted in finite time. Devel- 
opment of explosive instabilities is a characteristic feature of hydrodynamically unstable media (e.g. Pfirsch 
& Sudan 1993; Trubnikov et al. 1996). Plasma systems, like solar flares and TOKAMAKs disruption, often 
show explosive behavior (e.g. Cowley & Artun 1997). In these cases amplitude of perturbations A tends to 
infinity, A <~ (t — t) a , on Alfven time scale. 
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To illustrate explosive current disruption that may occur in the polar regions of GRBs we consider 
hydromagnetic collapse of a current carrying pinch (see also Liberman 1999). For simplicity we assume that 
collapse is non-relativistic and that the plasma obeys the ideal MHD equations. Since the size of the current 
carrying region m (Eq 30) is much larger than the light cylinder radius r s , we neglect the stabilizing effects 
of the poloidal magnetic field. The core is assumed to be cylindrically collimated and there is no velocity 
shear in the plasma rest frame. The core plasma is hot with temperature T, non-force-free and described by 
adiabatic index T a . Under these assumptions the governing equations are (all the quantities are measured 
in the plasma rest frame) 

d t p + —dvj(Tnpu) = 0, 

w 

d t B^ + d^(uB lp ) = 0, 

d t T + ud m T + (7 - l)—d„(wu) = 0, 
w 

p(d t u + ud^u)+d zo P+-^d^(wB^,= (El) 

where u is the (cylindrical) radial velocity. 

We will look for self-similar solutions of the system (El) where all the variables depend on the self-similar 
coordinate 

(E2) 

where R c (t) is the radius of the pinch at time t. Introducing scale parameter 

a(t) = Ml, (E3) 

where R c , > ro is the initial radius of the core, we parameterize the variables 

u = R c £U(0, p = p a 2x A(Z), Bfi = B^fia^H^), T = T a- 2X O(0- (E4) 
Assuming a special form of the velocity profile {/(£) = 1 we find 

£ = -l, M =-l, A = 7 -l (E5) 
with the following equations for the functions A, H p hi and G: 

H 4> { H 4>£)t 



e 2 A 
(QA)^ 
^A 



= Ci 

C 2 (E6) 



where Ci and C 2 are constants which can be set to unity. The radial dynamics of the pinch is then governed 
by 

where (i — 4ttT p /Bq is the plasma pressure parameter and time is measured in Alfvcn times, ta = 
B < f > fi/2 s /npoRo. The first term on the rhs of Eq. (E7) is responsible for the pressure support, while the 
second term is due to the magnetic field pinching. 
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In the spirit of our electro-magnetic approach we assume that the pressure forces cannot halt a collapse 
and set (3 = 0. Choosing constant density, A = const, we find that 

oc i (E8) 

The radial dynamics of the pinch is then governed by 

«"(*) = -J- (E9) 
a{t) 

Since the second derivative is negative, the scale factor a will become zero in finite time. This can be seen 
by direct implicit integration of eq. (E9) 

1 = \f^ Er f(V^h) (eio) 

where Erf is the error function. For a — > this gives 

, t= (Ell) 

Thus, the radius of a pinch becomes zero at t = Alfven times. 

Hydrodynamical description of the pinch collapse is applicable only for initial stages of the collapse. As 
the pinch contracts its radius will become smaller than ro- At this point the MHD approximation becomes 
invalid: there is not enough particles in the core to support the required current. As a result large inductive 
electric field E z > will develop. This will lead to particle acceleration and production of synchrotron 
radiation. 



F. Resistive collimation 

As we have discussed in Section 5.5 late collimation of relativistic outflows is hard to achieve due to 
effective relativistic suppression of lateral dynamics. In this section we show that collimation of relativistic 
outflows may be more efficient if we allow for dissipation of magnetic energy close to the axis and if the 
dissipated energy can decouple from the plasma flow, e.g. as radiation. If initially dissipation occurs near 
the axis, e.g. to enhanced resistivity brought about by current concentration, the resulting loss of magnetic 
pressure, which resists the inflow of plasma towards the axis, will propagate as a rarefaction wave away from 
the axis. After the propagation of the rarefaction waves plasma acquires lateral velocity which will leads to 
the pile-up of magnetic field near the axis and to faster radial expansion (toothpaste tube effect). 

As a simple model problem consider the structure of a core part of the flow in the flow rest frame. For 
didactic purposes we assume that the central core is cylindrically collimated and that there is no velocity 
shear in the plasma rest frame. In this case, the electric field vanishes and there is only the magnetic field 
due to a line current I': b = 21' /w' (w is the cylindrical radius in the plasma rest frame). 

Consider cylindrical line current RFF pinch in which at moment t = resistivity is turned on. Since 
there are no distributed charges, no distributed currents, equations of resistive RFF then become 



d t B = -curlE 
<9tE = curlB 



(Fl) 
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Eq. (Fl) then becomes an equation for cylindrical waves emitted by the line current. 

Mg^_^ = 27(t) *M 

w w 
For a given dependence I(t) its solution can be found by standard Green function method: 



-I 



„. „ (F3) 



Analytical solution exits, for example, for I{t) = Io<d(—t) + Iq(1 — t/ra)Q(t) where To is typical decay time 
for the axial current. We find 



27 (lnro-- ^ '-L-L j (F4) 



for t > r, and A z = 2I lnr for t < r. The fields are then given 




TO 



27, 

T() 

(see Fig. 11). 



E z = — In ((t - Vt 2 - w 2 )/m) (F5) 

Tn V / 



We can introduce a radial velocity with which electromagnetic energy is advected towards the resistor 
on the axis 

E z wln((t-VW^^)/w) 
Pro = -5- = , (i 1 D j 

Eq. (F6) has important implications. For \Jt 2 — w 2 ~ t electromagnetic velocity becomes of the order of 
the velocity of light and formally becomes larger than the speed of light (see Fig. 11). This illustrates a 
property of relativistic force-free plasma: there is no dynamical constraint that would preclude electric fields 
becoming larger than magnetic. Of course, as E — > B the applicability of RFF will break down and one 
need to use full relativistic MHD equations. One possibility then is that the plasma inertia will always keep 
the plasma from developing unphysical solutions. Alternatively, RFF dynamics will lead to "wave breaking" 
- creation of dissipative regions where particles will be accelerates (c/. Lyutikov & Blackman 2000). 
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Fig. 1. — Overall view of the model. 
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Fig. 2. — Current flow in the electromagnetic bubble. Current flow mostly along the axis, on the surface of 
the magnetic shell, along equator and close-up at the trailing part of the shell. Magnetic shell is preceded 
by the forward shock, typically r/T 2 ahead of it. Non-sphericity of the shell, which is of the order ~ 1/r 2 , 
is enhanced. 
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Fig. 3. — Magnetic "phase separation" near the central source. The currents supporting the strong magnetic 
field flow in the matter dominated phase inside the neutron star (or in the accretion disk in case of BH-torus 
system). The outside corona is magnetically dominated (c/. the Sun). 
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Fig. 4. — Evolution of strongly magnetized optically thick flows. Four-velocities of flows are given as 
functions of r/rj for L/M = 100 and different values of the parameter Flows start at r — rf with 
(3 = f3f; supersonic flows first accelerate as /?7 ~ r, reaching a terminal value given by the larger root of eq. 
(13), while subsonic decelerate initially as (3 ~ r~ 2 reaching asymptotic value given by the smaller root of 
eq. (13). 
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Fig. 5. — Magnetization parameter a. For supersonic flows (lower branch) the magnetization remains 
constant, reaching (7^ = (l — as r — > oo. Subsonic flows become strongly magnetized as they expand 

(upper branch); the magnetization parameter increases a ~ r 2 / 3 . 
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Fig. 7. — Development of pinching instabilities leads to core contraction, filamcntation, disruption of currents 
and particle acceleration by DC electric field. 
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Fig. 8. — Lorentz factor (dashed lines) and location of the contact discontinuity (solid lines) of the B&M 
solutions as a function of m for two choices of parameters: k = and k = 2. At the limiting value 
rrimax = k + 1 (point explosion cases) the contact discontinuity moves to \ — 00 • At the extreme case, for 
m = — 1, the contact discontinuity merges with the forward shock at \ = 1. 
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Fig. 9.- Lorentz factor (long dashed lines) and pressure (short dashed lines) as a function of \ f° r two 
choices of parameters: {to = 2,k = 0} and {m = 0,A; = 2}. Locations of the contact discontinuities are 
denoted by solid vertical lines. 



Fig. 10. — Variability of emission from relativistically moving emitters. A shell is moving with relativistic 
Lorentz factor T. In addition, primary emitters have "thermal" spread with a typical Lorentz factor < r. 
Each isotropic emitter produces a pulse of width A9 ~ l/(2r7r) when observed in the lab frame, while only 
emitters located within the angle 2/T may be seen by an observer. 




Fig. 11. — (a) Evolution of the toroidal magnetic field in the cylindrically collimated core due to dissipation 
of axial current (Eq. F5). At t = the axial current becomes dissipative with linearly decreasing total 
current, tq = 5, t = 0,1,3,5,6. Electromagnetic rarefaction wave propagates away from the axis, (b) 
Radial electromagnetic velocity for t = 1,3,4,4.5. When (3 r exceeds unity (in absolute value) force-free 
approximation breaks down. 



